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This paper presents precise large deviation estimates for solutions 
to stochastic fixed point equations of the type V = f(V), where 
f(v) = Av + g(v) for a random function g(v) — o(v) a.s. as v — > co. 
Specifically, we provide an explicit characterization of the pair (C, £) 
in the tail estimate P (V > u) ~ Cu~^ as u — > oo, and also present a 
Lundberg-type upper bound of the form P (V > u) < C(u)u~^. To 
this end, we introduce a novel dual change of measure on a random 
time interval and analyze the path properties, using nonlinear renewal 
theory, of the Markov chain resulting from the forward iteration of 
the given stochastic fixed point equation. In the process, we establish 
several new results in the realm of nonlinear renewal theory for these 
processes. As a consequence of our techniques, we also establish a 
new characterization of the extremal index. Finally, we provide some 
extensions of our methods to Markov-driven sequences. 

1. Introduction. 

1.1. Motivation and statement of the problem. Stochastic fixed point 
equations (SFPE) arise in several areas of contemporary science and have 
been the focus of much study in applied probability, finance, analysis of al- 
gorithms, page ranking in personalized web search, risk theory and actuarial 
mathematics. A general stochastic fixed point equation has the form 



where / is a random function satisfying certain regularity conditions and is 
independent of V. When f{v) = Av + B, where (A, B) is independent of V 
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and E [log A] < 0, this problem has long history starting with the works of 
Solomon (1972), Kesten (1973), Vervaat (1979), and Letac (1986). 

Tail estimates for solutions of general SFPEs have been developed by 
Goldie (1991) using implicit renewal theory, extending the fundamental work 
of Kesten (1973). Under appropriate regularity conditions, Goldie estab- 
lished that 

(1.2) P (V > u) ~ Cu-t as u -»• oo, 

where here and in the following, g{u) ~ h(u) means limu^oo g(u)/h{u) = 1. 
The constant C appearing in (1.2) is defined in terms of the tails of V 
and f(V), rendering it unsuitable for numerical purposes and for statistical 
inference concerning various parameters of scientific interest that involve 
the tails of V. Indeed, except in some special cases, the formulae for C 
presented in Goldie (1991) do not simplify to useful expressions. This issue 
is somewhat folklore and has been observed by several researchers including 
Yakin and Pollak (1998), Siegmund (2001), and Enriquez et al. (2009). Yakin 
and Pollak's work was motivated by likelihood ratio testing and change point 
problems in statistics, while the paper of Enriquez et al. was motivated by 
probabilistic considerations. 

The primary objective of this paper is to present a general alternative 
probabilistic approach to deriving the above tail estimates, yielding a char- 
acterization of this constant which — beyond its theoretical interest — is also 
amenable to statistical inference, and Monte Carlo estimation in particular. 
In the process, we draw an interesting connection between the constant C 
in (1.2) and the backward iterates of an SFPE described in Letac (1986). 
Specifically, we show that in many cases, this constant may be obtained 
via iterations of a corresponding perpetuity sequence, which has a compar- 
atively simple form and is computable. As this representation is explicit, it 
also resolves questions about the positivity of the constant which have been 
often raised in the literature and addressed via ad hoc methods. In addition, 
based on a slight variant of our method, we develop a sharp upper bound 
for the the tails of V, which has a form analogous to the classical Lundberg 
inequality from actuarial mathematics. 

The key idea to our approach is the observation that the distribution of 
the solution to the SFPE is frequently the stationary limit of a positively 
recurrent Markov chain and, hence, tail probabilities can be studied via ex- 
cursions within a regeneration cycle of this chain. After a large deviation 
change of measure argument, these excursions can be viewed as perturba- 
tions of a multiplicative random walk. While the multiplicative random walk 
itself can be handled via classical renewal-theoretic techniques, analysis of 
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perturbations requires the development of new methods in nonlinear renewal 
theory. 

The methods of this paper are quite general and apply to a wide class of 
SFPEs, including the following SFPEs which will be analyzed in some detail 
below: 



In the above, V is a real- valued random variable and independent of the 
other random variables on the right-hand side of these equations, while the 
random variables in the pair (A, B) or triplet (A, B, D) are allowed to be 
dependent on one another. Under the assumption that E [log A] < 0, a 
unifying feature of these equations is that they may all be written in the 
rough form V = AV+ "error," where the error term becomes negligible in 
the asymptotic limit for large V . 

It is also of considerable interest to study the case where E [log A] > 0. 
Borrowing the terminology from autoregressive processes, one may regard 
the case where E [log A] > as the explosive case while the case where 
E [log A] < as the stationary case. In a sequel to this work, we establish 
central limit theorems and rates of convergence results under appropriate 
scalings for the explosive case. Indeed, we establish that under some specific 
conditions, (V n /p n ) — > W as n — > oo with probability 1 (w.p.l) for a 
certain constant p. This limit random variable W will satisfy an SFPE, and 
the results of this article can be invoked to provide sharp tail estimates 
for W. Finally, the case where E [log A] = 0, referred to as the critical 
case in the literature, has been studied in detail in Babillot et al. (1997) 
and Buraczewski (2007) under fairly stringent conditions on the random 
variables B and D. However, methods of the current paper can be adopted 
to reduce some of these conditions. Hence, combining the results of this paper 
with known, but modified, results for the critical case and our results for the 
explosive case, one can obtain a complete description of the tail behavior 
for a wide class of random variables satisfying an SFPE encompassing (1.3)- 



We now turn to a more detailed description of specific examples that 
originally motivated this work. 



(1.3) 
(1.4) 
(1.5) 



V = (AV + B) + , 

V = Amax(D,V) + B, 

V = AV + BV 1/2 + D. 



(1.5). 



1.2. The ruin problem with stochastic investments. A fundamental prob- 
lem in actuarial mathematics — first introduced by Lundberg (1903) and 
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Cramer (1930) — considers the probability of eventual ruin for an insurance 
company. In this problem, the company is assumed to receive premiums at 
a constant rate c and incur claim losses which arise according to a Poisson 
rate. Thus, the capital growth of the insurance company is modeled as 

N t 

(1.6) X t = ct-Y,C, 

i=l 

where {N t } is a Poisson process and the loss sequence {£«} is assumed to 
be independent and identically distributed (i.i.d.). The ruin problem is con- 
cerned with the probability, starting from a positive initial capital of u, that 
the company ever incurs ruin. Since Xq = 0, this is equivalent to 

(1.7) := P (X t < -u, for some t > 0) . 
Cramer's classical result states 

(1.8) ^*(u)~C*e- ? " as u^oo, 

for certain positive constants C* and £, while the well-known Lundberg 
inequality states 

(1.9) **(«)< e - * u forallu>0. 

The constant £ is called the Lundberg exponent and the constant C* is 
explicitly characterized in (2.17) below. 

A natural criticism of Cramer's model is that the company's assets are 
assumed to be invested in a riskless bond yielding zero interest, while in 
practice these assets would typically be invested in a portfolio of stocks and 
bonds earning positive interest. One method to address this shortcoming 
is to introduce a stochastic investment process. To this end, let R n denote 
the returns on the investments in the nth discrete time interval, and let L n 
denote the losses from the insurance business during the same time interval; 
that is, 

L n = X n — X n -i, 

where {X t }t>o is as given in (1.6). The total capital Y n of the company at 
time n then satisfies the recursion 

(1.10) Y n := i? n y n _i - L n , n = l,2,..., Y = u. 

The ruin probability is now given by ^(u) := P (Y n < 0, some n), and by 
an elementary argument it can be seen that 



(1.11) 



ty(u) = P (C n > u, some n) , 
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where 

(1.12) C n := R^Lx + R l l R 2 l L 2 + ■ ■ ■ + (i?i ■ ■ ■ i?n) _1 L n . 

Using that R 2 ~ 1 L2 + • • • + (-R2 ■ ■ ■ Rn+i)~ 1 L n +i = C n , it follows from (1.12) 
that 

(1.13) C n = B + AC„-i where A= ^- and B = 

ill 

Hence, setting (sup ngZ £ n ) V = C, we obtain the SFPE 

(1.14) C = (AC + B) + , 

which is (1.3). Appealing to Goldie (1991)'s implicit renewal theorem yields 

(1.15) *(«) ~ Cu~^ as u -)• 00, 

obtaining a polynomial decay for ruin, which should be contrasted with the 
exponential decay in Cramer's original estimate. 

Variants on Goldie's result have been extensively studied in the actuar- 
ial literature; see, e.g., Paulsen (2008) and the references therein. For an 
extension to Harris recurrent Markov chains in general state space, see Col- 
lamore (2009). In spite of this extensive literature, we are not aware of a 
canonical method for characterizing the constant C in this and a variety 
of similar applications, and our paper fills this void. In practice, statistical 
estimates for C are of equal importance as those for £, and (1.15) is not 
very useful without a method for estimating these constants. Furthermore, 
as a by-product of our approach, we also obtain an analog of the classical 
Lundberg inequality for this problem, namely, 

(1.16) *(u) < C(u)u-Z for all u > 0, 

where C(u) converges to a positive constant as u — > 00. 

There is an interesting connection of this problem to the extremes of the 
well-known ARCH(l) and GARCH(1,1) financial models, which satisfy a 

related SFPE, namely V = AV + B. We will pursue this connection in more 
detail in Section 3, where we also apply the methods of this paper to yield 
a new, explicit characterization for the much-studied extremal index. 
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1.3. Related recursions and Letac's Model E. The models described in 
the previous section are special cases of a more general model introduced 
in Letac (1986). Specifically, Letac's Model E is obtained via the forward 
recursion 

(1.17) V n = A n max(D„, V n -i) + B n , n=l,2,... and V = v, 

where the driving sequence {(A ni B n , D n ) : n = 1, 2, . . .} is i.i.d., and we will 
assume in this discussion and throughout the paper that {A n } is positive- 
valued and E [log A] < 0. Iterating the above equation leads to the expression 



(1.18) V n = max [ II A i> V 

j=0 j=i+l k=l 



Y,Bi n A j + D k ]\A 3 

i=k j=i+l j=k 



with Bq = v. The sequence V n is a Markov chain, and one can show that 
under moment and regularity conditions, V := ]iai n -^ 00 V n exists almost 
surely (a.s.). A more compact expression is obtained in the form of an SFPE 
upon observing that (1.17) implies 

(1.19) V = Amsx(D,V) +B, 

which is (1.4). As we will discuss in more detail in the next section, there 
is also a corresponding backward iteration. Setting B* = Bi for i = 1, . . . , n 
and = v, then iteration of the backward equation yields that 



n+l 



(1.20) Z n = max V 



fc=i 



Y j B i \[A j + D k \[A j 



In fact, it can be shown that both the forward and the backward recursions 
satisfy the same SFPE, namely (1.19). 

Letac's Model E is quite general, encompassing a wide variety of recursions 
appearing in the literature. For example, setting D = — B/A in (1.19) yields 
V = (AV + B) + , which leads us back to the SFPE arising in our discussion 
of the ruin problem with investments. Similarly, the recursions arising in the 
ARCH(l) and GARCH(1,1) processes are easily seen to be special cases of 
(1.19) with D = and B strictly positive. Furthermore, we can also obtain 
the so-called perpetuity sequences arising in life insurance mathematics. Let 
B n be nonnegative, Vq = 0, and D n = for all n; then the backward 
recursion (1.20) simplifies to 



(1.21) Z n = B 1 + A X B 2 + A ± A 2 B 3 + ■ ■ • + (Ai • • • A n _i)S n , n = 1, 2, . 
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In a life insurance context, B n would represent future liabilities to the com- 
pany, and the quantity lim^oo Z n would then describe the discounted future 
liabilities faced by the company. The random variables A n are the "discount 
factors" which, as in the ruin problem with investments, may be viewed as 
the inverse returns of an appropriate financial process. 

Letac's Model E also includes reflected random walk, which arises in clas- 
sical queuing theory with the study of the steady-state waiting times for a 
GI/G/1 queue. Indeed, setting B = and D = 1 in (1.19), we obtain an 
SFPE corresponding to the forward equation 



which is seen to correspond with reflection at 1 of the multiplicative ran- 
dom walk nr=i Ai or, equivalently, to reflection at of the random walk 
S n := Y2i=i^°S^i- Using classical duality, the exceedances of {V n } can be 
related to maxima of the unreflected random walk. Similar arguments lead 
to further connections with classical risk models, autoregressive processes 
with random coefficients, and related problems. 

The rest of the paper is organized as follows. Section 2 contains notation, 
assumptions, background on Letac's principle and connections with Markov 
chain theory, and the main theorems of the paper. Section 3 is devoted 
examples and consequences of our main theorems. In Section 4 we describe 
our results in nonlinear renewal theory. These results serve as the essential 
technical tools that are needed for the proofs of the main results of this 
paper and seem to be of independent interest. Section 5 is devoted to proofs 
of the main theorems stated in Sections 3, while Section 6 contains the proofs 
of results on nonlinear renewal theory. An extension of Letac's Model E to 
nonlinear models is described in Section 7, while Section 8 is devoted to 
extensions for the Markov case, and Section 9 to some concluding remarks. 

2. Statement of results. 

2.1. Letac's principle and background from Markov chain theory. In the 
previous section, we presented forward and backward equations associated 
with the SFPE (1.19). In particular, the forward equation is generated by 
converting the SFPE (1.19) into the recursion (1.17), while the backward 
equation proceeds in similar fashion, but using an iteration backward in 
time. 

More precisely, suppose that we begin with the general SFPE 



n 



n 



(1.22) 




j=i i=j 



(2.1) 



V = f(V), 
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and assume that this equation can be written as a deterministic functional 
equation of the form 

(2.2) V = F(V,Y), 

where F : R x R d -)• R is assumed to be deterministic, measurable, and 
continuous in its first component. For convenience of notation, we will often 
write F y (v) = F(v,y), where y £ R d . For example, the SFPE (1.19) may 

be written as V = F Y (V), where F(v,Y) = Amax(D,v) + B and Y = 
(A, B, D). 

Let v be an element of the range of F, and let {Y n } be an i.i.d. sequence 

v 

of random variables such that Yi = Y for all i. Then the forward equation 
generated by the SFPE (2.2) is defined by 

(2.3) V n (v) = F Yn oF Yn _ 1 o...oF Yl (v), n = l,2,..., V = v, 
while the backward equation generated by this SFPE is defined by 

(2.4) Z n (v) = F Yl oF Y2 o...oF Yn (v), n = l,2,..., Z = v. 

Note that the backward recursion need not be Markovian, but for every 
v and every integer n, V n (v) and Z n (v) are identically distributed. Letac's 
(1986) principle states that when the backward recursion converges to Z a.s. 
and is independent of its starting value v, then the stationary distribution 
of V n is the same as the distribution of Z. This is precisely described in the 
following lemma. 

Lemma 2.1. Let F : R x R d — > R be deterministic, measurable, and 
continuous in its first component; that is, F Y is continuous, where Y : R d — > 
R for some d € Z + . Suppose that lim n ^. 00 Z n (f) := Z exists a.s. and is 
independent of its initial state v. Then Z is the unique solution to the SFPE 
(2.2). Furthermore, V n {v) =^ V as n — > oo, where V is independent of v, 
and the law of V is same as the law of Z . 

Starting from the general SFPE (2.2), our basic approach will be to gen- 
erate the forward recursive sequence (2.3) and study the evolution of this 
process. In the next section, we will impose additional assumptions that 
will ensure that the forward recursive sequence viewed as a Markov chain 
possesses a regenerative structure, which we now describe. 

For the remainder of this section, let {V n } denote a Markov chain on a 
general state space (§, S), and assume that this chain is aperiodic, countably 
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generated, and irreducible with respect to its maximal irreducibility mea- 
sure, which we will denote by 92. Let P denote the transition kernel of {V n }. 
In Section 2.2 onward, we will take {V n } to be the Markov chain generated 
by a forward recursion, but for the remainder of this section we will assume 
it to be an arbitrary chain. 

Recall that if {V n } is aperiodic and c^-irreducible, then it satisfies a mi- 
norization, namely 

(M) 5l c {x)v{E) < P k (x, E), Vi6§,£G5, k€ Z+, 

for some set C with (p(C) > 0, some constant 5 > 0, and some probability 
measure v on (<S,§). The set C is called the "small set" or "C-set." It is no 
loss of generality to assume that this C-set is bounded, and we will assume 
throughout the article that this is the case. 

An important consequence of (A4), first established in Athreya and Ney 
(1978) and Nummelin (1978), is that it provides a regeneration structure for 
the chain. More precisely, these authors established the following. 

Lemma 2.2. Assume that (A4) holds with k = 1. Then there exists a 
sequence of random times, < Kq < K± < ■ ■ ■ , such that: 

(i) Kq, K\ — Kq, K2 — Ki, . . . are finite a.s. and mutually independent; 

(ii) the sequence {Ki — : i = 1, 2, . . .} is i.i.d.; 

(iii) the random blocks {Vk 1 _ 1 , • • • , V^-i} are independent, i = 0, 1, . . .; 

(iv) P (V Ki e E\d Kl -i) = v(E), for all EeS. 

Here and in the following, $ n denotes the a— field generated by Vo, V\, . . . , V n . 
Let Tj := Ki — denote the ith inter-regeneration time, and let r denote 

a typical regeneration time, i.e., r = t\. 

Remark 2.1. Regeneration can be related to the return times of {V n } 
to the C-set in {M) by introducing an augmented chain {(V n ,r] n )}, where 
{r] n } is an i.i.d. sequence of Bernoulli random variables, independent of {V n } 
with P (rj n = 1) = 5. Then based on the proof of Lemma 2.2, we may iden- 
tify Ki — 1 as the (i + l)th return time of the chain {(y n ,r] n )} to the set 
C x {1}. At the subsequent time, then has the distribution v given in 
(M), independent of the past history of the Markov chain. Equivalently, the 
chain {V n } regenerates with probability 5 upon each return to the C-set. In 
particular, in the special case that 5 = 1, in which case the chain {V n } is 
said to have an atom, the process {V n } regenerates upon every return to the 
C-set. 
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Note that if the Markov chain is also recurrent, then r < oo a.s. How- 
ever, in this article, we will also consider transient chains which may never 
regenerate. In that case, we will set r = oo. 

As consequence of the previous lemma, one obtains a representation for- 
mula (Nummelin (1984), p. 75) relating the stationary limit distribution of 
V to the behavior of the chain over a regeneration cycle. Set 

r-l 

(2.5) iV« = E 1 K-)(V;), 

n=0 

where Vq ~ v. Thus, using Lemma 2.2 (iv), we see that N u describes the 
number of visits of {V^} to the set (u, oo) over a typical regeneration cycle. 
The representation formula alluded to above is given in the following lemma. 

LEMMA 2.3. Assume that (A4) holds with k = 1. Then for any nGl, 

(2.6) P(F >„) = ^1. 

2.2. Results in the setting of Letac's Model E. We now turn to the core 
results of the paper. In this discussion, we will focus primarily on Letac's 
Model E, namely, 

(2.7) V = F Y (V), where F Y (v) = A max(£>, v) + B 

for Y = {A, B, D) G (0, oo) x M? (although we will later generalize the 
method). As described previously, this model is sufficiently general to incor- 
porate the main applications described in the introduction. Let {V n } and 
{Z n } denote the forward and backward recursions generated by this SFPE, 
respectively. In addition to these sequences, we also introduce two related 
backward recursive sequences, which will play an essential role in the sequel. 

The associated perpetuity sequence: Let {(A n , B n , D n )} ne i + be the i.i.d. 
driving sequence which generates the forward recursion and let Aq = 1 

and Bq ~ v for v given as in Lemma 2.2 (iv), where Bq is independent of the 
driving sequence {(A n , B n , D n )} ne % + . Now consider the backward recursion 

(2.8) ZV=F$ a...ofj(0), n = 0,l,..., 
where 



(2.9) 



== £ + | ■ 
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By an elementary inductive argument, it follows that 

n _ 

The sequence {Zn^} will be called the perpetuity sequence associated to {V n } 
and is the backward recursion generated by the SFPE (2.9). 

The conjugate sequence: Let Dq := —Bq and D* := —AiD-i — Bi for 
i = 1,2,..., and consider the backward recursion 

(2.11) ZV=F$o...of£(0), n = 0,l,..., 
where 

(2.12) F^\v):=jimn(D*,v) + ^. 
It follows by induction that 

The sequence {z'tf*} will be called the conjugate sequence to {V n } and is the 
backward recursion generated by the SFPE (2.12). 

Further notations: Let A denote the multiplicative factor appearing in 
(2.7), and define 

A(a) = logEL4 Q ] and A(a) = log A(a), Va. 

Let /x denote the distribution of Y := (log A, B, D), and let /x Q denote the 
a-shifted distribution with respect to the first variable; that is, 

(2.14) n a (E) := [ e ax dfi(x,y,z), \/E G B(M. 3 ), Va£l, 

Je 

where here and in the following, B(M. d ) denotes the Borel sets of M rf for any 
d G Z+. Let E a [•] denote expectation with respect to this a-shifted measure. 
For any random variable X, let £(X) denote the probability law of X. Also, 
write X ~ £(X) to denote that X has this probability law. Let supp(X) 
denote the support of X. Also, given an i.i.d. sequence {Xi}, we will often 
write X for a "generic" element of this sequence. 

For any function h, let h! denote its first derivative, and so on. Also let 
dom (h) denote denote the domain of h and supp (h) the support of h. Fi- 
nally, let 3n denote the cr-field generated by the Letac's forward sequence 
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Vq,Vi,. . . , V n , where the process {V n } is obtained from the recursion (2.7). 
We now state the main hypotheses of this paper. 

Hypotheses: 

(H±) The random variable A has an absolutely continuous component 
with respect to Lebesgue measure, and satisfies the equation 
A(£) = for some £ G (0, oo) n dom (A'). 

(H 2 ) E[|£|*] < oo and E[(A\D\)^] < oo. 

(H 3 ) P (A > 1, B > 0) > or P (A > 1, B > 0, D > 0) > 0. 

Next we turn to a characterization for the constant of decay. To describe 
this constant, we will compare the process {V n } generated by Letac's forward 
equation (2.3) to that of a random walk, namely 

5„:=log^i + --- + log^ n , n = l,2,..., 

where {log^4j} is an i.i.d. sequence with common distribution (ia- In Sparre- 
Andersen's generalization of the classical ruin problem, one studies ruin for 
the random walk {S n }, that is, 

*$>*(u) := P (S n < -u, for some n £ Z+) , 

and it is well known that, asymptotically, 

(2.15) ^*(n)~C*e- ?u as u -+ oo, 

where E[^] = 1. The constant C* may be described in various ways, but 
a characterization due to Iglehart (1972) is obtained by first setting 

(2.16) r* =inf{n> 1 : S n < 0} , 
and then defining 

1 - E[e ?s ^*] 

(2 - 17) c * = WOW*}' 

We will refer to C* as the Cramer-Lundberg constant. 

Our aim is to develop an analogous characterization for the constant of 
tail decay for the sequence {V^}. Expectedly, it will involve the return times 
r, where r is a typical regeneration time of the Markov chain in the sense of 
Lemma 2.2. In Lemma 5.1 below, we will show that {V n } is indeed a Markov 
chain satisfying appropriate path properties for this lemma to hold. 

Returning to the central result of the paper, recall that {Z ( ? ] } and {Z^} 
are the associated perpetuity and conjugate sequences, respectively, and let 
Z( p ) and denote their a.s. limits. In Lemma 5.5 below, we will show 
that these limits exist and that Z^ > Z^ a.s. 
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Theorem 2.1. Assume Letac's Model E, and suppose that (Hi), (H2), 
and (H3) are satisfied. Then 

(2.18) lim u f P (V > u) = C 

for a finite positive constant C which is given by 



(2-19) C = Eg 



1 

£V(OE[r] 



(z^-z^yi {T ^ } 



Remark 2.2. Comparing (2.19) with (2.17), we obtain alternatively 
that C = C*6E^[(Z^ - Z^)h {T=oo} ] , where 

It should be noted that the expectation in the previous theorem is com- 
puted in the ^—shifted measure, and for this reason, the event {r = 00} 
occurs with positive probability. Moreover, the constant C obtained in The- 
orem 2.1 can be decomposed into three parts, of which we would like to em- 
phasize the Cramer- Lundberg constant C* and Eg [(Z( p ) — Z( c ))^l j r=00 |] . In 
the proof, it will be seen that the constant C* results from the "large-time" 
behavior of the process which, under the rare excursion to the high level u, 
resembles the random walk S n = Yli=i 1°§ A- Indeed, when (B,D) = (0, 1), 
then the constant C reduces to the Cramer-Lundberg constant C*. In con- 
trast, the constant Eg[(z( p ) - Z^fl 

{ r =oo}J results from the "small time" 

behavior of the process. 

For the elementary recursion defined by f(v) = Av + B, Enriquez et al. 
(2009) also provide a probabilistic representation for C but their character- 
ization is complex and is not easily amenable for statistical computation. In 
contrast, we obtain a comparatively simple form for the constant, namely 
Eg [(z( p ) — Z^Y^-{t=od}] j although the exact form of this constant will de- 
pend on the particular recursion under study. However, in a wide variety 
of applications, the conjugate term Z^ will actually be zero on {r = 00}, 
and so the integrand involves only the perpetuity term, Z^ , which has the 
simple form of (2.10) and whose terms can be shown to converge to zero at 
a geometric rate. 

Returning to the ruin problem with investments described in Section 1.2, 
the forward recursive sequence is given by 

(2.20) V n = (A n V n -! + B n ) + , n = l,2,..., V = v. 
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First, to determine the regeneration time r for this process, begin by ob- 
serving that the minorization condition (M) holds with 

(2.21) 6=1, C = {0}, and v = £(B + ). 

In other words, {V n } has an atom at {0}. Consequently, by Remark 2.1, we 
may take r — 1 to be the first return time of the process {V n } to the origin, 
starting from an initial state Vq, where Vq ~ v = £(B + ). Next we turn to 
identifying the conjugate sequence. Upon setting D{ = —Bi/Ai in (2.13), we 
obtain 

(2.22) *' = -KAE^l 

\ n =l i=0 U V 

But the terms on the right-hand side all have the form Y^=o " " " A) 

for some n, which upon multiplication by (Aq ■ ■ ■ A n ) yields 

n 

(2.23) Y / B i {A i+1 ---A n ), n = 0,l,.... 

i=0 

Now for to be nonzero, we would need the quantity in (2.23) to be 

negative for some n. But (2.23) is just the nth iteration of the SFPE V = 
AV + B with forward recursive sequence {V„}, and from the definition of 
{V^} and {Vy, we have 

inf{n : V n < 0} = inf{n : V n < 0} := r - 1, 

since these two processes agree up until their first entry into (— oo,0], at 
which time the reflected process {V n } is then set equal to zero. Hence we 
conclude that, conditional on {r = oo}, = and so the constant C 
is determined entirely by the perpetuity term Z&\ and conditioning on 
{r = oo} simply removes those paths for which the process {V^} returns to 
the origin in its ^-shifted measure. 

Next we supplement the previous estimate with a Lundberg-type upper 
bound. First define the perpetuity sequence 

(2.24) Z^ = \V Q \ + ± (l ^ i +A j A|) l { ^ } - 

i=l % 

Theorem 2.2. Assume the conditions of the previous theorem. Then 



(2.25) 



P (V > u) < C(u)u-t, for all u > 0, 
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where, for certain positive constants C\{u) and C2(u), 

(2.26) C {u) = EM Z^] sup \ e^ z {zd{u) + C 2 (u))\ <oo for all u. 

Furthermore, 

1 cr 2 

(2.27) CAu) -)• — and C 2 (u) -)• H ^ as n ^ oo, 

m 777/ 

w/iere m := Eg [log A] and cr 2 := Varg (log A). 

The constants Ci(u) and C*2(n) will be identified explicitly in (6.43) below. 

2.3. Generalizations of Letac's Model E. While the results of the previ- 
ous section were limited to Letac's Model E, the method of proof is actually 
more general, applying to a wide class of recursions. For example, by a slight 
modification of the proof, it is also possible to analyze polynomial recursions 
of the form 

(2.28) V = B k V + B k _M k ~ 1)/k + ■■■ + BxV 1 ^ + B , 

where the random variable B k is assumed to satisfy the same conditions 
as the random variable A in Letac's model, and the vector (Bq, . . . , B k ) 
takes values in (0, co) k+1 . We obtain the same asymptotic results as before, 
although the form of the constant C is more complicated. We will return 
to this example and some further generalizations in Section 7, and to an 
extension to a Markovian driving sequence in Section 8. 

3. Examples and applications. In this section, we develop a few ex- 
amples and draw a connection to the related problem of estimating the 
extremal index. 

3.1. The ruin problem with stochastic investments. Returning to the ruin 
problem described in Section 1.2, we recall that in this case the SFPE spe- 
cializes to f(v) = (Av + B) + . 

In the context of this problem, it is interesting to point out that there is a 
close connection of our approach to methods based on duality. Specifically, 
to obtain the above SFPE, we began by studying the risk process in (1.10), 
which leads, via Letac's forward recursion, to the Markov process 



(3.1) 



V n = {A n V n -i + B n ) + , n = l,2,..., y = o, 
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which we will analyze in the proofs of our theorems below. Alternatively, 
we could have arrived at the process in (3.1) using a duality technique as 
described in Asmussen and Sigman (1996). In particular, the main result 
of that article states that the "content process" {V^} is dual to the "risk 
process" {Y n }. In particular, their risk process is obtained by inverting the 
function f(v; a, b) = (av + b) + for fixed a, b, and then studying the forward 
recursion generated by f*~, which is easily computed to be the function {Y n } 
in (1.10). Hence by Asmussen and Sigman (1996), Theorem 3.1, 



(3.2) 



P (V > u) = P (Y n < 0, for some n G Z+) := *(«). 



This shows — by an alternative approach — that the probability of ruin may 
be equated to the steady-state exceedance probability of the forward recur- 
sive sequence {V^}. 

Appealing now to Theorem 2.1 of this paper, we obtain (2.18) with 



(3.3) 



C = 



£A'(0E[r 



where r — 1 is the first return time of {V n } to the origin and Bq ~ £(B + ); 
cf. the discussion following (2.21). Note that the last expectation may be 
simplified for £ > 1, since Minkowskii's inequality followed by a change of 
measure argument yields the upper bound 



Bi 



2^ A n ... 4, 



An - ■ ■ A; ' {T=0 ° } 



< E 



\B\t 



oo 

i+x;p(r><-i) 



1=1 



Thus the constant C may be explicitly related to the return times of the 
process {V n } to the origin. Since {V n } is geometrically ergodic (Lemma 5.1 
below), these regeneration times have exponential moments, and so this last 
expression is also finite. 

Finally, if £ < 1, then an analogous result holds by applying the deter- 
ministic inequality 



(3.4) 



\x + y| < 



V(x,y) £l 2 , £G (0,1], 



in place of Minkowskii's inequality. For some specific examples of investment 
and insurance processes, see Collamore (2009), Section 3. 
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3.2. The ARCH (1) and GARCH(1,1) financial processes. From a math- 
ematical perspective, the ruin problem described in the previous example is 
closely related to a problem in financial time series modeling, namely, the 
characterization of the extremes of the G ARCH (1,1) and ARCH(l) financial 
processes. 

In particular, in the GARCH(1,1) model of Bollerslev (1986), the loga- 
rithmic returns on an asset are modeled as 

(3.5) R n = o- n ( n , n=l,2, 

where is an i.i.d. Gaussian sequence of random variables, and where 
{o- n } satisfies the recurrence equation 

(3.6) a 2 n = a + bo- 2 n _ x + cR 2 n _ x 

for certain positive constants a, b, and c. The quantity a n is called the 
stochastic volatility. Substituting (3.5) into (3.6) yields the recurrence equa- 
tion 

(3.7) V n = A n V n -i + B n , n = l,2,.... 

with V n = a 2 , A n = (6 + c^_ 1 ), and B n = a. We note that this same 
recurrence equation also arises in the well-known ARCH(l) model of Engle 
(1982). 

By comparing this SFPE with the SFPE of the previous example, we see 
that we are dealing with essentially the same process as before, except that in 
this case the random variables B n are always positive, meaning that we may 
no longer take r to be the first return time to {0}. To determine a correct 
choice of r we need to verify, as before, the minorization condition (M) of 
Section 2.1. To this end, note that in the GARCH(1,1) model, V(v) := Av + 
B has a density, which we call h v . Since h v (x) is monotonically decreasing 
for large x, 

inf h v (x) = ho(x) for sufficiently large x. 

ve[o,M] 

Hence 

0(E) := / I inf h v (x) ) dx 
Je \ve[o,M] J 

is not identically zero, and then P(x,E) > l[ ^(x)C'(E). Finally, upon 
normalizing v so that it is a probability measure, we obtain (M) with k = 1 
and 5 the normalizing constant. Thus, according to Remark 2.1, we see that 
regeneration occurs w.p. S upon return to the set C = [0, M\. 
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Once this modification has been made in the choice of r, we may repeat 
the calculation in (3.3) of the previous example to obtain an explicit char- 
acterization of the constant C in Theorem 2.1. As in the previous example, 
we obtain that the conjugate term is zero on {r = oo}, in this case 
because the process {V n } is always positive, and so we indeed obtain (3.3) 
but with a different r. 

3.3. The extremal index. While the previous examples and estimates 
have dealt with stationary tail probabilities, the methods can also be ap- 
plied to other related problems. The extremal index — which measures the 
tendency of a dependent process to cluster — is defined for a strictly station- 
ary process {V^} by the equation 

(3.8) P ( max V„ < u n ) = e~ et , 

\l<i<n J 

for any given t, where {u n } is chosen such that 

(3.9) lim n {l-P(> n >u n )} =t. 

The constant in (3.8) is called the extremal index. As argued in Rootzen 
(1988), p. 380, this quantity may be written in a Markovian setting as 

(3.10) 9= lim P(K>n,somen<r) 

where N u denotes the number of exceedances above level u occurring over 
a regeneration cycle. We will show in this article, using the methods of 
Theorem 2.1, that 

3.11 9 = r L . J , 

v ' E[r*] 

where r* = inf{n > 1 : S n < 0} for S n = Y27=i l°S^i- The above expression 
provides a considerable simplification of known theoretical estimates for G, 
as given, e.g., for the special case of the GARCH(1,1) model in de Haan 
et al. (1989) and, moreover, this last expression is seen to hold in the general 
setting of Letac's Model E. 

4. Some results from nonlinear renewal theory. The main tools 
needed in the proofs of Theorems 2.1 and 2.2 will involve ideas from nonlin- 
ear renewal theory, which we now present in detail. 
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Recall that the crux of the proof of Theorem 2.1 will be to utilize the 
representation formula of Lemma 2.2, namely, 

E [N u ] 



P (V > u) 



EM 



Thus we will need to estimate E [N u ] , the number of exceedances above 
level u which occur over a regeneration cycle. To study this quantity, we will 
employ a dual change of measure, defined as follows. First let (A\, B\,Di) ~ 
fi£, where fi£ is the ^-shifted measure defined previously in (2.14). If V\ > u, 
then let (A n , B n , D n ) ~ jjl for all n > 1. Otherwise let (A2, B2, D-i) ~ 
and then for V2 > u set (A n , B n , D n ) ~ \x for all n > 2, while for V2 < u 
set (A3, B3, D3) ~ fj,£, etc. In summary, the dual change of measure can be 
described as follows: 

(4.1) £ (lo g A,, B „,A,) = {« tlZk" T " 

where T u = inf{n : V n > u} and ^ is defined as in (2.14). Roughly speaking, 
this dual measure shifts the distribution of log A n on a path terminating at 
time T u , and reverts to the original measure thereafter. Let Ej> [•] denote 
expectation with respect to the dual measure described in (4.1). 
We now focus on two quantities: 

(i) the overjump distribution at the first time V n exceeds the level u, 
calculated in the ^-shifted measure; and 

(ii) the expected number of exceedances above level u which then occur up 
until regeneration, calculated in the original measure. 

Note that in the ^-shifted measure, E^ [log A] > and hence V n t 00 w.p.l 
as n — > 00 (c/. Lemma 5.2 below). Consequently, 

V n = A n max (D n , V n -i) + B n 

implies 

V n ~ A n V n -i for large n. 

In other words, for large u the process {V n } will ultimately resemble a pertur- 
bation of multiplicative random walk in an asymptotic sense. Consequently, 
the problem described in (i) may be viewed as a variant of results from 
nonlinear renewal theory (c/. Woodroofe (1982)). 

Lemma 4.1. Assume Letac's Model E, and suppose that (Hi), (H2), and 
(H3) are satisfied. Then in the dual measure (4.1), 



(4.2) lim Pt I —>y 

u->oo \ u 



T U <T 



) = P e (y > y) 
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for some random variable V. The distribution of this random variable V is 
independent of the initial distribution of Vq and is described as follows. If 
A 1 is a typical ladder height of the process S n = Y]™—] log Ai in the ^-shifted 
measure, then 

1 f°° 

(4.3) P 5 ( log V > y) = —j-j, / P ? (A 1 >z)dz, for all y > 0. 

E £ [ A \ Jy 

While Lemma 4.1 follows rather immediately from known results in non- 
linear renewal theory, this is not the case for the problem stated in (ii) above. 
Here we would like to determine E53 [N u |3t„]j *-e., the expected number of 
exceedances, in the dual measure, which occur over a typical regeneration 
cycle prior to the regeneration time r. We provide a precise estimate for 
this quantity next. In the following result, note that we start the process at 
a level vu, where v > 1, and so the dual measure actually agrees with the 
original measure in this case. 

Theorem 4.1. Assume Letac's Model E, and suppose that (Hi), (H 2 ), 
and (H3) are satisfied. Then for any v > 1, 



(4.4) lim E 



V 

= V 

u 



= U(\ogv), 



where U(z) := ^ neN fJ*^(— 00, z) and ha is the marginal distribution of 
-log A. 

Roughly speaking, the function U may be interpreted as the renewal func- 
tion of the random walk — S n = — Y17=i ^S-^i- More precisely, if the dis- 
tribution of log A is continuous, so that the open interval (—00, z) could be 
replaced with the closed interval (—00,2:], then U(z) agrees with the stan- 
dard definition of the renewal function of {— S n }. 

We now supplement the estimate in the previous theorem with an upper 
bound. 

Theorem 4.2. Assume that the conditions of the previous theorem. Then 
there exist finite positive constants C\(u) and C2(u) such that 

(4.5) E S [N u |ff T „A(r-i)] < (c^log^) +C 2 (u)j l {Tu<r} , all u, 

where C\(u) and C^u) are the positive finite constants that converge as 
u — > 00 to to" 1 and 1 + (a 2 /m 2 ), respectively, where m := E^ [log A] and 
a 2 ■= Var^ (log A). 
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We emphasize that this last theorem is crucial for obtaining the Lundberg 
upper bound in Theorem 2.2. 

Finally, using the previous results, it is now possible to state an extension 
of Lemma 4.1, which will be particularly useful in the sequel. First let 

Q u :=Es, [N u \d Tu A( T -i)] 
In the following, we write Pz,(A) for Ej> [1^] . 

Theorem 4.3. Assume the conditions of the previous theorem. Then 
conditional on {T u < t}, 

(4.6) Q u =► U(]ogV)V~* as u -»• oo, 

where V is given as in (4.2) and (4.3). That is, 

lim ( Q u < y\ T u <t) = ([/(log V)V^ < y) , for all y > 0. 
5. Proofs of the main results. 

5.1. Five preliminary lemmas. We begin by presenting five preparatory 
lemmas, which will be needed in the proofs of the main theorems to be 
proved in Section 5.2. The first two lemmas establish path properties of the 
Markov chain {V^}, obtained via forward iteration of the SFPE. Recall that 
a Markov chain satisfies a drift condition if 

(V) / h(y)P(x, dy) < ph{x) + /31 c (x), for some p e (0, 1), 

where h is a function taking values in [1, oo), (3 is a positive constant, and 
C is a Borel subset of R. In the next theorem, we work with the assump- 
tion that the process {V n } is nondegenerate, namely that we do not have 
P (Fy(v) = v ) = 1. This assumption is subsumed by (#3), but we state it 
here separate hypothesis. 

Lemma 5.1. Assume Letac's Model E, and let {V n } denote the forward 
recursive sequence corresponding to this SFPE. Assume that (Hi) and (H2) 
are satisfied and that the process {V n } is nondegenerate. Then: 

(i) {Vn} satisfies the drift condition (V) with C = [—M,M]. 

(ii) {V n } is (p -irreducible, where ip is the stationary distribution of {V n }. 

(iii) {V n } satisfies the minorization condition (M) with k = 1 under both 
the measure p and the measure p^. Furthermore, the set [—M,M] is petite 
for any M > 0. 



u 



-{T u <t}- 
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(iv) {V n } is geometrically ergodic. Moreover E [e €T ] < oo, where r is the 
inter-regeneration time of the process {V n } under any minorization (A4) 
where the C-set is bounded. 

Proof, (i) Using the defining equation V n = A n max (D n , V n -i) + B n , 
we obtain the inequality 

(5.1) \V n \ < A n \V n ^\ + (A n \D n \ + \B n \) . 

Now hypothesis (Hi) implies that E [A a ] < 1 for all a G (0,£). Fix a G 
(£ A 1) and Vq = v. Then from the deterministic inequality (3.4), we obtain 

(5.2) < E[A a ]v a + E[\B 1 \ a ] +E[(A 1 \D 1 \) a ]. 

Thus (D) holds with h(x) = \x\ a + l, p = (E [A a ] + 1) /2, and C = [— M, M], 
where M is chosen sufficiently large such that 

^ (1 - E L4 Q ]) M Q > E[|Bi| a ] + E[(Ai|I>i|) a ] . 

(ii) To verify that {V^} is irreducible, it is sufficient to show that if <p(A) > 
0, then U(x, A) > for all x, where 

oo 

U(x,A) :=Y,P n {^A) 
n=i 

(cf. Meyn and Tweedie (1993), Proposition 4.2.1). Now to identify 99, note 
under our hypotheses that the conditions of Proposition 6.1 of Goldie (1991) 
are satisfied. Hence the backward recursion (1.20) converges and is indepen- 
dent of v. Thus, by Lemma 2.1, it follows that {V^} is stationary and the 
law of the limiting random variable V is the same as the law of Z. Now 
choose ip = £(Z). Since {V n } is aperiodic, it follows by Theorem 1.3.1 and 
(1.16) of Meyn and Tweedie (1993)) that 

(5.3) sup \P n (x,A) - ip(A)\ -> asn^oo, Vx G R 

(where we have used (i) to verify their condition (1.14) with h(x) = |x| a + l). 
Hence 



(5.4) U(x,A) > lim P n (x,A) = <p(A) > 

n— >oo 

for any ^-positive set A. 
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(iii) We note that the main issue is to verify that the minorization holds 
with k = 1 (rather than for general k). 

We will show that for any v G R, there exists an e- neighborhood Q3 e (t>) 
such that 



for some positive constant 5 and some probability measure v, both of which 
will typically depend on v. Since some such interval *B e (v) will necessarily 
be 93-positive with (p = £(V) as in (ii), this will imply the existence of a 
minorization with k = 1. 

Set v* = inf {v : P (D < v) = 1} € (—00, 00]. We will consider three dif- 
ferent cases, namely v < v*, v > v*, and v = v* . 

If v < v* , then ~P (D > v) > 0. Since P (-D > w) is nonincreasing as a 
function of w, P (D > v + e) > for some e > 0. For this choice of v and e, 
let *B e (f) be the required e- neighborhood in (5.5). Now if Vq G ^8 e (w) and 
Di > v + e, then max(Z?i, Vb) = D\ and hence Vi = A\D\ + Thus, for 
any initial state Vb G 2$e(v), 



Note that the first term on the right-hand side is independent of Vq. Taking 
v to be the probability law of (A\D\ + B\) conditional on {D\ > v + e} and 
5 = P (D > v + e), we obtain (5.5). 

Next suppose v > v* . Then for some e > 0, P (D > v — e) = 0, and 
for this choice of v and e, let 5S e (w) be the e-neighborhood in (5.5). Then 
P (V\ = A\Vq + B\ I Vb G 33 e (-u) ) = 1. Hence, to obtain a minorization for V\ 
in this case, it is sufficient to derive a minorization for V\{v) := A\v + B\ a.s. 
To this end, begin by observing that since A\ has an absolutely continuous 
component, so does the pair (A±, A\v + B\). Let ^'"^ denote the probability 
law of {A\,A\v + B\). Then by the Lebesgue decomposition theorem, 7^) 
can be decomposed into its absolutely continuous and singular components 
relative to Lebesgue measure, and the continuous component satisfies 



where I denotes Lebesgue measure on R 2 . Hence by construction, there exists 



(5.5) 



P(v, E) > 8l^ e{v) v{E), for all E G B(R), 



(5.6) 



Vi = (A i D 1 + Bi) l{ Dl>1)+e } + Vil{ Dl <i, +e }. 




a rectangle where 



dl 



is bounded away from zero; that is, 



(5.7) 



(x,y) > S > 0, for all (x, y) G [a, b] x [c,d]. 



dl 
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Now suppose w G and let denote the probability law of V\{w) := 

A\w + B\. Then V\(w) — V\(v) = A\(w — v), and hence (5.7) yields 

—j^{x,y)>S, for all (x,y) E [a,b] x [c + w*,d — w*], 

where w* = b\w — v\. Notice that the constant b was obtained by observing 
that the first component describes the distribution of A, where supp(^4) C 
(0, 00), implying that < a < b < 00. Since w* ], as w — > v , it follows that 
for sufficiently small e, there exists a subinterval [d , d'] of [c, d] such that 

(5.8) — ^ — (x, y) > 5, for all (x,y) G [a, b] x [c',d'] and all u> G ® e (v). 
Hence the minorization (5.5) holds with 

v(E)= f ( inf ^L(x,y) j d/(x,y), for all E 1 G £>(M), 
JrxE \™e!Be(v) at / 

and by (5.7), this measure is not identically equal to zero. 

It remains to consider the case where v = v* . First observe that, similar to 
(5.6), P (Vi = AiVq + Bi \V G %>e(v*), D 1 < v* - e) = 1. Hence, we clearly 
have 

(5.9) V 1 = (A 1 V + Si) l {Dl <^- e} + Vil {Dl> ,*_ £} . 
Now define 7( e ' 1 '*) on any Borel set E C R 2 to be 

7 M*)(£) = p((Ai,Aiv* + Bi) 6£,Di <v*-e) ^ ^ V *\E) as e -> 0, 

where the last step follows from the definition of j^ v *\ since P (D± < v* — e) t 
1 as e — >■ 0. Hence, since 7^*) has an absolutely continuous component, there 
exists an e > such that ^ e < v *^ also has an absolutely continuous compo- 
nent. Now apply the previous argument with 7^) replaced with ^ e < v *^ for 
this choice of e to obtain the corresponding minorization for this case. Thus 
we have obtained (5.5) for all the three cases. Note that the above compu- 
tations hold regardless of whether we are in the original measure \x or in the 
^-shifted measure //£. 

Finally, to show that [—M,M] is a petite set, note that [—M,M] C 
U»€[-mm] %$e(v)( v )i where *B e ( v )(v) is a small set and hence is petite. Thus 
there exists a finite subcover of petite sets, and then by Proposition 5.5.5 of 
Meyn and Tweedie (1993), [-M,M] is petite. 



TAIL ESTIMATES FOR HOMOGENEOUS NONLINEAR RECURSIONS 



25 



(iv) Since [— M, M] is petite for any M > 0, it follows from (i) and Meyn 
and Tweedie (1993), Theorem 15.0.1, that the process {V n } is geometrically 
ergodic. It remains to show that, regardless of (5, C,u) in the minorization 
(Ai), E [e eT ] < oo, where r is the inter-regeneration time under (A4). Now 
since C is bounded, it follows by Theorem 15.2.6 of Meyn and Tweedie (1993) 
that C is /i-geometrically regular with h(x) = \x\ a + 1. Consequently, letting 
K denote the first return time of {V n } to C, we have 

(5.10) F(t) := supE [t K \V = v] < oo 
for some t > 1. 

Consider the split-chain (see Nummelin (1984), Section 4.4 and Remark 
2.1 above). Originating from the C-set, this chain has the 1-step transition 
measure v(dy) w.p. 5 and the 1-step transition kernel P(x, y) := (P(x, dy) — 
v(dy))/(l — 5) w.p. (1 — 5). Thus, conditional on either one of these transition 
kernels, it follows from (5.10) that T u (t) < oo and Tp(t) < oo, where T u 
denotes conditioning on V\ ~ u, and Tp(t) := sup^g^ E [t K \ Vq = v, Vi ~ P] . 
Hence by Remark 2.1, 

(5.11) E [f\ < T v (t) [d + 5(1 - 5) n (T P (t))^ , 

where 5 is the constant appearing in (M). Finally observe by a dominated 
convergence argument that Tp(t) J, 1 as t \. 1, and thus for sufficiently small 
t > 1, rp(i) < (1 — <5) _1 . For this choice of t, the series on the right-hand 
side of (5.11) is convergent, as required. □ 

In the next result, we establish a critical result concerning the transience 
of the process {V n } in its ^-shifted measure. 

Lemma 5.2. Assume Letac's Model E, and let {V n } denote the forward 
recursive sequence corresponding to this SFPE. Assume that (H\), (H2) and 
(H3) are satisfied. Then under the measure 

(5.12) V n /■ +00 w.p. 1 as n — > 00. 

Thus, in particular, the Markov chain {V n } is transient. 

Proof. We begin by showing that the chain {V n } is transient in the 
^-shifted measure. To this end, first note by (#3) that the set [M, 00) is 
attainable with positive probability, for any positive constant M and any 
initial state. Thus <p([M, 00)) > for all M. Consequently by Meyn and 
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Tweedie (1993), Theorem 8.3.6, it is sufficient to verify that for some positive 
constant M, 

(5.13) P ? (V n < M, for somen € Z + \V > 2M) < 1. 
To establish (5.13), first note by definition of {V n } that 

(5.14) V n > A n V n -i - \B n \, foralln, 
and iterating this equation yields 

(5.15) - Vn A >V - W n , where W n := ^ + ■ ■ ■ + lD " 1 



Ai---A n ~ u A 1 A x ---A n 

We will now show that {W„} converges in distribution to a proper non- 
degenerate random variable. To this end, note that since {(A n , B n )} ne i + is 
an i.i.d. sequence, W n has the same distribution as 

n \r>\ n 
i=l j=i+l 

where {V n } is a Markov chain. Moreover, {V^} satisfies hypotheses (Hi) and 
(H2) and is nondegenerate. Consequently by Lemma 5.1 (iv), {V n } converges 
to a proper distribution and hence so does {Wn}- 

Next observe that [log A] = A'(£) > 0; thus, the i.i.d. sequence {log A n } 
has a positive drift in the £— shifted measure. Consequently, setting S n — 
Yli=i ^°sAi we obtain P (S n < 0, for some n) < 1. Therefore 

(5.16) P f (Ai ■ • • A n > 1, for all n) = p > 0. 

Now assume that (5.13) fails. Then by (5.15) and (5.16), we would have 

(5.17) P f (W n > M, for some n) > p, for all M. 

But {W n } is nondecreasing and converges to a proper distribution. Hence 
(5.17) is impossible, so we conclude (5.13). 

Next observe by Lemma 5.1 that the set [—M,M] is petite for every 
M > 0. Hence [— M, M] is uniformly transient (Meyn and Tweedie (1993), 
Theorem 8.0.1) and, consequently, the expected number of visits of {V n } to 
[-M, M] is finite for all M. It follows that \V n \ t +00 w.p. 1 as n ->■ 00. But 
we cannot have V n 1 —00 with positive probability since, by definition, we 
have as a lower bound that V n > A n D n + B n for each n, where {(B n , D n )} 
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forms an i.i.d. sequence of random variables. Consequently V n f +00 w.p.l 
as n — > 00. □ 

Prior to stating the next lemma, recall that T u := inf {n : V n > u}, and 
that Ej) [•] denotes expectation with respect to the dual measure described 
in (4.1). We now compare an expectation in the original measure to this 
expecation in the dual measure, following a typical regeneration cycle from 
its initial state with Vb ~ v, where v is the measure described in Lemma 2.2 
(iv), and terminating at the next regeneration time r. 

Lemma 5.3. Assume the conditions of Lemma 5.2. Let g : M°° — >■ R be 
a deterministic function, and let g n denote its projection onto the first n + 1 
coordinates; that is, g n (xo, X\, • • •) = g(xo, . . . , x n ). Then 

(5.18) E[0 T _i(K o ,...,V r _i)]=E a) [g T ^(Vo,...,V T ^)e-^l {Tu<T} 

+ Es, [g T -i(V , . . . , F T _i)e-^l {Tu > r} 

Proof. If 

£ (log A k ,B k ,D k ) = 

then it can be shown using induction that 

(5.19) E [g n (V , V n )] = E € [g n (V , V n )e^ s " 

Now (5.18) follows by conditioning on {T u = m,r = n} and summing over 
all possible values of m and n. □ 

Next we establish a critical result which links Letac's forward iteration of 
the SFPE to the backward equations. 

Lemma 5.4. Assume Letac's Model E, and let {V n } denote the forward 
recursive sequence corresponding to this SFPE. Let {Z^} and {Z„^} denote 
the associated perpetuity sequence and the conjugate sequence, respectively, 
assumed to have the initial values described above in (2.8) and (2.11), and 
set Aq = 1. Then for any n G N, 



fi^ for k = 1, . . . , n, 
fi for k > n, 



(5.20) (z<f)-Z( c ))l { ^ >0} = Z„l { ^ >0} , where Z n :-- 



Vn 

Aq ■ ■ ■ A n ' 
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Proof. It follows by an inductive argument that for all n G N, 



n n 



(5.21) K = max ( II Aj, \J 

=0 j=i+l k=l 



n n 



Y, >'>■ II Aj + D^Aj 



i=k j=i+l 



j=k 



where -Bo : = Vq. Hence 
(5.22) 
where 
(5.23) 

Now 
(5.24) 



V„. 



E 



k=i 



fe-i 

E 

.i=0 



Bi 



min (0, 9Jt n ) , 



A ---Ai A ---A k . 



4 P) =^ ) o-o^(0) = Et A T- 



=0 



A ---Ai 



(c) 

To obtain a similar expression for Z n , observe again by induction that 

k 



(5.25) 4 C ) =min J] 



A 



An ■ ■ ■ A 

u 1 k=0 



E 



+ 



Dt 



=o 



A) • ■ ■ ^ Aj • • • A fc 



Then substituting .Dq = — B and D* = — AA — -Bi, « = 1,2, ... , yields 
that with Tl n given as in (5.23), 

(5.26) Z^ =min(^),0,93T n ) . 

Comparing these expressions with (5.22), we conclude that 

(5.27) Z n = -mm(0,Tl n ). 

Finally observe that if Z n > 0, then it follows from the previous equation 
that < max (z ( n v) , Z n p) - 9JT n ) and hence 



Z n l {2nm = max (0, ZM, z£> - 9JT n ) = Z^-Z^ 



c) 



□ 



In the next lemma, we study the convergence of the sequence {Z n } and 
relate it to the random variable Z^ defined just prior to the statement of 
Theorem 2.2. 
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Lemma 5.5. Assume the conditions of Lemma 5.2, and let {Z n } be de- 
fined as in (5.20). Then in ^-measure, {Z n } has the following regularity 
properties: 

(i) Z n — > Z a.s. as n — > oo, where Z is a proper random variable supported 
on (0, oo). As a consequence, Z^> and Z^ are well defined and Z^ > Z^ 
a.s. 

(ii) < oo. Moreover, for all n and u, 



(5.28) 



{n<r} 



< Z^ 



and 



< Z^\ 



Proof. We begin by establishing (ii). Since 

(5.29) \V n \ < A n \V n ^\ + (A n \D n \ + \B n \) , 

the process {|V^|} is bounded from above by {R n }, where Ro = \Vq\ and, 
for each n £ Z+, 

(5.30) Rn = B n + A»-Rn-i, where B n = {\B n \ + A n \D n \) . 
Iterating the previous equation yields 



(5.31) 



R n = Y,Bi [] Ai, n = 0,l,..., 

i=0 i=j+l 



where Bq = | Vb | - Since \V n \ < R n for all n, we may now apply (5.20) to 
obtain that {|Z n |} is bounded from above by the perpetuity sequence 



(5.32) 

and hence 
(5.33) 



W n := (A 1 ---A n )~ 1 R n = J2 



Bi 



i=0 



Ay -AC 



\z, 



Bi 



n|l{n< r } < W n l{ n<T y < ^ - _*^. 1 {r>i} : ~ Z{P \ 



and similarly for Zt u 1{t u < t }- 

It remains to show that E[(Z( p ))^] < oo. To this end, observe that if 
£ > 1, then by Minkowskii's inequality followed by a change of measure 
argument, 



(5.34) 



[Z^f] ^ < E [#] 1/? ( 1 + £ P (t > i ~ l) Vt ) < 



oo, 



i=l 
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where finiteness is obtained from (H2) and Lemma 5.1 (iv). On the other 
hand, if £ < 1, then an analogous result is obtained using the deterministic 
inequality (3.4) in place of Minkowskii's inequality, which then completes 
the proof of (ii). 

To establish (i), set Z_ = liminf „^oo Z n and observe that if P(Z_ < 0) = 
p > 0, then by Fatou's lemma 

(5.35) p < liminf E [l fi <( J = liminf P (V n < 0) , 



where, in the last step, we have used Lemma 5.4 to replace Z n by V n . But 
by Lemma 5.2, V n f +00 w.p.l as n — > 00, and thus (5.35) is impossible. We 
conclude Z_ G (0, 00). 

Next observe by definition of {V n } that if V n -i > 0, 

(5.36) A n V n -! - \B n \ <V n < A n Vn-i + {\B n \ + A n \D n \) . 

Let S(n) denote the indicator function on the event {V/ > 0, for all I > n}. 
Then by Lemma 5.2, hin^^oo 5{n) = 1 a.s. By iterating the left- and right- 
hand sides of (5.36), we obtain that for any m > n, 

m n 

(5.37) \V m -(A n+1 ---A m )V n \S(n)< £ ^ II A, 

i=n+l i=j+l 

where Bi is defined as in (5.30). Hence for any m > n, 

00 £ 

(5.38) Z m -Z n 5(n)< ^ - ^ . 

i=n+l ^ * 

Now recall that [log ^4] = A'(£) > 0. Hence the process 

S n := - log Ax log A n + log B n 

has a negative drift, and using (Hi) we have E^ [(A -1 )"] < 00 for some 
a > 0. Hence by Cramer's large deviation theorem (Dembo and Zeitouni 
(1998), Section 2.2), it follows that 

(5.39) P 5 (S n > n (-E c [log A] + e)) < e~ tn for all n, 
for some constants e G (0, Eg [log A]) and t > 0. Thus 
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where a := exp {— E^ [log A] + e} G (0, 1). Then by the Borel-Cantelli lemma, 

< 541 » ^(i^ >0 " i0 ) =a 

Since lim ra _ >00 S(n) = 1 a.s., it follows as a consequence of (5.38) and (5.41) 
that Z n converges a.s. to a proper random variable on (0, oo), which we 
denote by Z. 

Since the above argument also holds for the recursion f(v) = Av + B, it 
follows that Zn^ converges a.s. to a random variable taking values in (0, oo). 
Then by Lemma 5.4, the limit in the definition of Z^ must also exist and 
we must have that - Z^ = Z > a.s. □ 



5.2. Proofs of the main theorems. The key to these proofs is the following 
result, which shows that k^E [N u ] behaves asymptotically as a product of 
two terms, the first being influenced mainly by the "short-time" behavior 
of the process, while the second term is determined by the "large-time" 
behavior. We will later identify this second term as a classical limit related 
to the random walk process S n = Y17=i 1°S^- 

PROPOSITION 5.1. Assume Letac's Model E, and suppose that (Hi), 
(H2), and (H3) are satisfied. Then 



(5.42) lim n^E [N u ] = E € 



&1 



{r=oo} 



lim Ejj) 

u— >co 



u 



T U <T 



PROOF. Since N u = g(Vo, . . . , V T -\) for some function g, it follows by 
Lemma 5.3 that 

(5.43) E [N u ] = E s [N u e~^ . 
Moreover by the definition of Z n in (5.25), 

(5.44) e-* 5T «l { T„<r} = ({^■■■A Tu )^vll {Tu<T} )v-t 

= (4.1{T.<r } ) V^. 

Combining the last two equations yields 

(5.45) E [N u ] = Ej, [(4 u l {Tu<r} ) N U V~^ 
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Consequently, by conditioning on $t u a(t-i) an d rearranging terms, we have 
that 



(5.46) 
where 

(5.47) 

Note 
(5.48) Ed 



u«E [JV u ] = E 5) [(4 it l {Tu<T} ) Q u 



Q u := Ex, [N u \$t u a(t- 



i). 



u 



1 {T u <t}- 



^tJ-{t u <t} 



Qu 



— Ed 
+ Ed 



ZiQul{ n <T u <T} 
(^T„ 1 {T u <r} - ^nl{n<T„<r}) Qu 



Now take the limit first as u — > 00 and then as n — > 00. 

We begin by analyzing the first term on the right-hand side. By Theorem 
4.2 and by the definition of Q u in (5.47), {Q u } is uniformly bounded in u. 
Also, by Lemma 5.5 (ii), E^ [|Z n |?l| ri<Tu<r }] < 00 for all n. Consequently 
for any given n, 



(5.49) 



Ed 



nQu\l{n<T u <T} < KE^ \Z n \ l{ n < Tu<T } 



<K' 



for finite constants K and K'. In the middle term, we have replaced Ed [•] 
with E^ [•], since these will be the same on the set {n < T u < r}. Hence an 
application of the dominated convergence theorem yields 



(5.50) 



lim Ed 

u— >oo 



Zi,Qu^-{n<T u <T} 

= Ed 
= Ed 



%ri lim 1 {n<T u <r}^ [Qu |3n, T u < t] 

Zi lim l {n<Tu<r} j Ed [Q] , 



where the last step of (5.50) follows by Theorem 4.3 (and, in particular, the 
independence of the limiting distribution Q on the initial distribution of the 
process {14}, which is here taken to be the distribution of V n for some fixed 
n). In the second expectation on the right-hand side, we have exchanged the 
limit and expectation, which is justified based on the weak convergence of 
Qu to Q and the uniform boundedness of {Q u } (obtained from Theorem 4.2 
and the definition of Q u ). 

To identify the first expectation on the right-hand side of (5.50), observe 
once again that Z n l{ n < Tu<T j is the same in the ^-shifted measure as it is in 
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the dual measure, since the dual measure agrees with the ^-shifted measure 
up until time T u . Hence 



E 



Zi lim 1 



{n<T u <r} 



Zi lim 1 



{n<T u <r} 



E, 



zii 



{r=oo} 



since an elementary argument yields that T u f oo a.s. as « — > oo. 
Substituting into (5.50) and now taking the limit as n — > oo yields 



lim lim E 



Z i,QvX{n<T u <r} 



= lim E^ 



{t=oo} 



[Q] 



= E^l {r=oo} ]E s [Q], 



where the last step follows from the dominated convergence theorem and 
Lemma 5.5 (i) and (ii). Finally, observe by Theorem 4.3 that 

Ej, [Q] = lim Ej, [Q u ] (P (T u < t))- 1 



= lim Ej) 

M— S>00 



lim Ej) 



E S [N u \d Tu A(r-l)] 



1 {T u <t} 



(p (r„<r)) 



-1 



u 



T U <T 



[In the second equality, we have used that N u = on the set {T u > r}.] 
Substituting the previous equation into (5.50) yields 



(5.51) 



lim lim Ej> 

n— »oo u— >co 



Z iQul{n<T u A(T-l)} 



E ? [Zl {r=oo} ] lim Ej, 



T u <t 



Going back to (5.46) and (5.48), we see that the proof of the proposition 
will now be complete, provided that we can show that 



(5.52) 



lim lim E35 



Z tJ-{T u <t} 



Z ll{n<T u <T}) Q 



0. 



But recall once again that {Q u } is uniformly bounded. Thus, to establish 
(5.52), it is sufficient to show 



(5.53) 

and 

(5.54) 



lim lim Et 

n— >oo u— >oo 



lim lim Ec 

n— >co u— >oo 



i Z T u 1 {Tu<r} ~ Z&l{ n <T u <T}) 







= 0. 
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Note that in these last expectations, we have again replaced Ej> [•] with 
[■], since these expectations involve random variables on {T u < r}, and 

on that set these expectations are actually the same. 

To establish (5.53), first apply Lemma 5.5 (ii). Namely observe that the 

integrand on the left-hand side is dominated by 2(Z^)^, which is integrable. 

Hence, applying the dominated convergence theorem twice, first with respect 

to the limit in u and then with respect to the limit in n, we obtain 



(5.55) E 4 



(lim Z% \ {Tu<r} - lim Z|l {T=oo} ) 



0, 



where, in the last equality, we have used that T u f 00 as u — > 00. This 
establishes (5.53). The proof of (5.54) is analogous, so we omit the details. □ 

Next we identify second term on the right-hand side of (5.42) by relating 
it to the classical Cramer-Lundberg constant. First recall that r* is a typical 
return time of the random walk S n = Ya=i ^-i t° the origin, while r is a 
typical regeneration for the process {V n } under study. 

Lemma 5.6. Assume Letac's Model E, and suppose that (Hi), (H2), and 
(H3) are satisfied. Then 



5.56 I J lim E s 



N, 



Vt u 



u 



T U <T 



= C*, 



where C* is the Cramer-Lundberg constant defined in (2.17). 

Proof. Define 

W n = (log A n + W n ^) + , for n = 0, 1, . . . . 

Then {W n } is a random walk reflected at the origin, and since E [log -4] < 0, 
it is well known that this process is a recurrent Markov chain which converges 
to a random variable W whose distribution is the stationary distribution. 
Also, by Iglehart (1972), Lemma 1, 

(5.57) lim u-tp (W > logn) = C* , 

where W := lim^oo W n . 

Set r* = inf {n € Z + : W n = 0}. Since {M^n} has an atom at the origin, 
t* + 1 is a regeneration time of the Markov chain {H^}. Hence by the 
representation formula in Lemma 2.3, 

e[at*1 

(5.58) P(W > logu) = L ul 



E[r*] ' 
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where 

T* 
71=1 

Let ha denote the marginal distribution of log A, and set 

HA*(E)= f e^ x dfi A (x), V£eS(l). 
Je 

Let T* = inf {n 6 Z + : W n > logn}, and define the dual measure 

With a slight abuse of notation, let E33 [•] denote expectation with respect 
to this dual measure. Then by a change of measure, analogous to Lemma 
5.3, we obtain 



(5.60) E[N*]=u-tEs> 
= u _$ E s 

which is equal to 



N* e -t.(W T *-\ogu) 
E S [K |^ AT .] e-^- logu )l {T ,< T , } 



'J) 



-€(W T .-logu) 



T* < T * 



P S (T*<r*). 



To complete the proof, notice that the multiplicate random walk process 
{exp(W n )} also satisfies the conditions of Theorem 4.3, and the random 
variable inside the last expectation converges to the same weak limit as the 
corresponding quantity for the process {V n }. (Note that we write "T* < r*" 
while we write "T u < r," since r is a regeneration time for the original 
process, while r* + 1 is a regeneration time for the reflected random walk 
process.) Thus by Theorem 4.3, ( Ji(u)/ J2{u)) — > 1 as u — > oo, where 



u 



T u < T 



Ji(u):=Ea, Ev[N u \$ Tu 
J 2 (u) := E s [Ej, [N* |ff T . AT ,] e -^ T£ -io g u) 
Consequently we conclude 



(5.61) lim E s 



V T 
u 



lim 



u^oo Ps, (T* < t*) 
C*E [r*] 



lim — ^ by (5.57) and (5.58). 

u^oo P^, (T* < T*) y V ' V ' 
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Finally observe that lim u _>oo (T* < r*) = P33 (T* = oo), and by a change 
of measure argument, 

oo 

(5.62) Pj, (r* = oo) = 1 - P ? (r* < oo) = 1 - Eg [l {r * =n} ] 

?1=1 

OO 

= l-^E[ e ^l {r , =n} ]=l-E[e^]. 



n=l 



The required result then follows from (5.61) and (5.62). 



□ 



PROOF of Theorem 2.1. By Lemma 2.3, P (F > u) = E[A^]/E[r], 
and hence by Proposition 5.1, 



(5.63) lirn u^P (V > u) = Eg Z^l {r=oo} 

■ (E[r]) _1 lim E s 



u 



T U <T 



By Lemmas 5.4 and 5.5 (i), Z = — Z^ a.s., and by Lemma 5.6, the 
last limit on the right-hand side of (5.63) may be identified as 

C*E [r*] 1 



1-E[e^*] ^'(0 



by (2.17). 



□ 



PROOF of Theorem 2.2. A repetition of (5.46) and (5.47) yields 

'V T . x 



(5.64) u^B [N u ] = E 



'T 

u 



< E 



(Z<P))*(Ci(«)log 



+ C 2 («) 



u 



where the last step follows from Theorem 4.2 and Lemma 5.5 (ii). Moreover, 
since Vr u > u, we obviously have 



(5.65) 



Ci(u) log 



u 



+ C 2 {u) 

<sup(e-^ (zC 1 (u) + C 2 (u))\, 



which is bounded. Substituting this deterministic bound into the right-hand 
side of (5.64) establishes the theorem. The limiting values of C\(u) and 
C2{u) are obtained by a further application of Theorem 4.2. □ 
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Returning to the extremal index described in (3.11), recall that 



(5.66) 6= lim 

u— >oo 



P {V n > u, for some n < r) 



lim 

E [N u ] E [N u ] 

Moreover one can show, similar to the proof of Proposition 5.1, that 

Vt. x 



E [1 



{Tu<t}\ 



lim ««E [l {Tu<T} ] =E 5 Z*l {T=oo} 



lim E35 

U— >00 



L {r«<r} 



-i u 

U 



T U <T 



To identify the last term on the right-hand side, let T* and r* represent 
the corresponding terms in the reflected random walk model {W n } (as in 
the proof of Lemma 5.6). Then by applying Lemma 4.1 to the process {V n } 
and to the multiplicative random walk {exp(W / ra )}, we see that Vt u /u and 
log Wt* /logn converge to the same weak limit. Hence 



(5.67) lim Ex, 

u — ^oo 



U 



T u <t 



lim Ej, 

M— >• oo 



-£,{W T *-\ogu) 



T* < T* 



But the quantity on the right-hand side is just the Cramer-Lundberg 
constant C*. To see that this is the case, note that if T* j£ r*, then the 
unconstained random walk S n = Y17=i^°S^-i returns to the set (— oo,0], 
and starting from this new level, the overjump distribution will have the 
same limiting distribution as it had when starting its first cycle from the 
origin. Thus, using the fact that the dual and ^-shifted measures are actually 
the same in this case, we obtain 



lim Et 



-«W T .-]ogu) 



T* < T* 

^ u — ' 



lim E^ 



^(^-logu) 



T* < oo 



lim u^P(T* < oo) = C*, 



where the second equality follows from a change of measure argument and 
the fact that P^(T* < oo) = 1. Thus we conclude that 



(5.68) 



lim^E [l {Tu<T} ] = C*E^ \ztl {T=oo} 



Now by Proposition 5.1 and Lemma 5.6, we also have that 



(5.69) 



lim ti ? E [N u ] 



E, 



{t=oo} 



C*E [r*] 
l_ E [e^r*] 



Substituting (5.68) and (5.69) into (5.66), we conclude (3.11). 
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6. Proofs of the results from nonlinear renewal theory. First we 
turn to the proof of Lemma 4.1, which extends a classical result from nonlin- 
ear renewal theory (given in Woodroofe (1982), Theorem 4.2) to the setting 
of our problem. 

Proof of Lemma 4.1. By definition of Z n in (5.20), 

(6.1) V n = (A 1 ... A n ) (Z n l {2n>0} + Z n l feo} ) • 

Now by Lemma 5.2, V n t oo w.p.l under the measure fi^. Thus T u < oo a.s., 
and at this exceedance time, we obviously have l{z n <o} = ®- Thus Vt u = V^ u 
for 

(6.2) K = (A 1 ---A n )Y n , n = l,2,..., 

where Y n = Z n on {Z n > 0} and Y n = 1 otherwise, and thus Y n is everywhere 
positive. We begin by showing that 

(6.3) -± V as u -> oo, 

u 

where has the distribution described in (4.3). 
To this end, note by (6.2) that 

(6.4) logV^ = S n + 5 n , 

where S n := J27=i^°S^i an d S n := logy n . Hence {logl/^} may be viewed 
as a perturbed random walk where, under jUg-measure, {S n } has a positive 
drift and the sequence {S n } has the property that {(log Ai, Si) : i = 1, . . . , n} 
is independent of log Aj for all j > n. This puts us in the setting of classical 
nonlinear renewal theory. 

Next observe by Lemma 5.5 and \,g <Q y — > a.s. that 5 n := logl^ 
converges a.s. to a proper random variable, and hence S n /n — > a.s. as 
n — y oo. Thus {S n } is slowly changing. By Woodroofe (1982), Theorem 4.2, 
it follows that (6.3) holds, and hence this equation also holds with Vp u in 
place of Vt u ■ 

To show that the result holds conditional on {T u < r}, let Vq ~ u, where 
v is given in Lemma 2.2 (iv); let Kq = 0; and let K\, K2, ■ ■ ■ denote the 
successive regeneration times. For each i, let Ri denote the distribution of 
(Vt u /u) conditional on the event that T u 6 [Ki, Ki+i). By independence 
of the regeneration cycles, {Ri} is an i.i.d. sequence of random variables. 
Consequently, the conditional distribution of Vt u given that T u G [Ki, i^j+i) 
is the same for all i. Thus, the conditional distribution of Vt u /u given {T u < 
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t} must agree with the unconditional distribution of Vr u /u, and the proof 
is complete. □ 

Next we turn to the proofs of Theorems 4.1-4.3. An important part of 
the proofs will be the study the process {V n } as it returns from the rare 
set (u, oo). To do so, we will introduce a new barrier at a lower level it*, 
where t G (0,1), and divide the trajectory into two parts, namely that 
occuring prior to the event V n G [—if*, u% and that occuring after this time. 
Intuitively, the process {V n } will closely resemble a multiplicative random 
walk away from the set [—it*, it*], and so a critical aspect of the proofs will 
be to characterize the behavior of {V n } after it returns to [—it*, it*], but 
prior to regeneration. Indeed, it is precisely in this region that {V n } fails to 
resemble the random walk process. In the next proposition, we show that 
the behavior of the process in this critical region may be neglected in an 
appropriate sense, and we also provide a reasonably sharp estimate for the 
number of visits above the level it which arise after returning to [— n*, it*]. 

For any v > 0, first define 

K(v) = inf{n : \V n \ < v}. 

PROPOSITION 6.1. Assume Letac's Model E, suppose that {Hi), (H 2 ), 
and (H3) are satisfied, and assume that Vo > u. Then for any t G (0,1), 
there exist finite positive constants a, M, and p G (0, 1) such that 



(6.5) E 



T-l 

Yl 1 {Vn>u} 

n=K(u t ) 



< A(u). 



The constant A(u) is characterized as follows. Set V\ := A\ max (Di, Vq) + 
|J3i|. Then 

«-«(!-*) f _ ) 

(6.6) A(n) := — - { 1 + u~ at sup E w [r] E M [Vf] } < 00 

1 _ P { we[-M,M] J 

(where E w [•] denotes expectation conditional on Vq = v). 

Note that we have dropped the dependence on the dual measure in this 
proposition, since we assume that Vo > u, and hence the entire trajectory 
will take place in the original measure. 

PROOF of Proposition 6.1. Recall by Lemma 5.1 (ii) that {V n } satisfies 
a drift condition; namely, for some a > 0, 

(6.7) EhvX \Vn-i = v] < p\v\ a + P1 Cm (v), 
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where p G (0, 1), Cm '■= [-M, M], and M and (3 are constants. 

We will divide the proof into three steps. In the first step, we will study 
the number of exceedances above level u which occur in an excursion begin- 
ning at time K (u*) (i.e., when the process {|Vn|} first falls below the level u* 
where t < 1) and ending at time K(M) (i.e., when this process first enters 
the set Cm '■= [—M,M]. In the next step, we then consider the number of 
exceedances above level u which occur between time K(M) and the actual 
regeneration time. Combining these two estimates in Step 3, we obtain the 
desired upper bound. 



Step 1: For any t G (0, 1), 
(6.8) 



E 



K(M) 
n=K(u t ) 



< 



u -a(l-t) 
1-P ' 



Proof. Suppose Vq = v, where \v\ > M. Then iterating (6.7) yields 
(6.9) E [\V n \ a l {K(M)> n } \Vo = v]< p n \v\ a , n = 0, 1, . . . , 

and hence 



(6.10) E 
Consequently, 



t{V n >u}t{K{M)>n} 



Vq = v 



< u- a p n \v ia 



(6.11) 



E 



'K(M)-l 

1 {Vn>u} 

n=0 



V n = v 



<u~ a \v\ a {l + p + p 2 + ■■■) 



Then by the previous equation and the strong Markov property, 



(6.12) E 



'K(M)-l 

Y 1 {V n >u}\VK{u' ■ 

n=K{u>-) 



<u- a \V K(ut) \ a (l + p + p 2 + ■■■), 



Since iVR-^td < u l by definition, this establishes (6.8). 
Step 2: We have 

r-l 



(6.13) 



E 



Yl 1 {Vn>u} 

n=K(M) 



U 



< < sup K,,. [r] E M [K 

we[—M,M] 



1-p 
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Proof: It suffices to show that for any v G [— M, M], 

(6.14) E [N U \V =v]< ^— \ sup E w [r] E M [V?] \ . 

1 ~ P [u>e[-M,M] J 

To this end, introduce the augmented chain described in Remark 2.1. 
Namely consider the process {(V^,,?7 n )}, where {r/ n } is an i.i.d. sequence of 
Bernoulli random variables with P {rj n = 1) = 5 and P (rj n = 0) = 1 — 5. 
As discussed in Remark 2.1, we may then identify r as the return time of 
this augmented chain to the set € := C x {1}, where C and 5 are obtained 
from the minorization condition (M) (which holds by Lemma 5.1 (hi)). We 
assume without loss of generality that C C Cm, and define Wl = Cm x {0, 1}. 

With a slight abuse of notation, let P denote the transition kernel of the 
augmented chain {{V n ,rj n )}. Then introduce the taboo transition kernel 

G P(x,E) := f l G c(y)P(x,dy), 
Je 

where G and E are Borel subsets of 1 x {0, 1}. Then by the last-exit de- 
composition (cf. Meyn and Tweedie (1993), Section 8.2), 

(6.15) £ P"(v, E) = OT P"(v, E) + Y J f c P fc (v, dw) a „P"- fc (w, E). 

Now set E = (u, oo) x {0, 1} and fix v = (v, q), where v < u and q £ {0, 1}. 
Then sum (6.15) over all n £ Z+. On the left-hand side of (6.15), the term 
£ _p n (v, (u, oo) x {0, 1}) describes the probability that regeneration does not 
occur during the first n time increments and that V n G (u, oo). Thus 

oo 

(6.16) ^eP"(v>,oc)x{0,l}) 

n=l 

oo 

= J^E [l {v . n>u} l {r> „ } |y = v] = E [N U \V = v] . 

n=\ 

Next consider the second term on the right-hand side of (6.15). Applying 
Tonelli's theorem to interchange the order of summation and integration, 
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and then interchanging the order of summation, we obtain 

oo n— 1 „ 

( 6 - 17 ) EE/ ^ fc (v,dw) OT P™- fc (w,(n,oo)x{0,l}) 

„ oo / oo 

= / E^ fe ( v ' dw ) E OT ^ fc ( w >K^)x{o,i}) 

7 « fc=l \n=k+l 



-f I 



V € P k (v,dw) m P(w,dz) 
/we(OT-c) Jzean c fe=1 



]T aKP n (z,(n,oo) x{0,l}) 



Vn=0 



For the last term in parentheses, note that gtj{P n (z, (u, oo) x {0, 1}) describes 
the probability that {V n } avoids the set [—M,M] during the first n time 
increments and that V n £ (it, oo). Hence setting z = (z,r) yields, by (6. If), 



oo oo 

(6.18) ^ OT P n (z,(n,oo)x{0,f})=E £ 

1 {V n >u} 1 {K(M)>n} 



n=0 


_n=0 






"A"(M)-1 




= E 


E 1 {Vn>«} 


V = z 




ra=0 





V n = z 



< 



u- a \z\ a 



1 - 



Substituting this inequality into the previous equation, we see that the right- 
hand side of (6.17) is bounded above by 

(6.19) ^E/ eP k (v,dw)([ m P(w,dz)\z\A. 

Note that with w = (w, s), 

[ m P(w,dz)\z\ a <B[V 1 a \V = M] , 

where V\ := A\ max (D\, Vq) + \P>\\. Moreover, for the remaining integral in 
(6.19), we have by definition that 

/ c P k (v, dw) = P(r>k,V k £ [-M, M]\V = v) . 
Jm-t 

Substituting the last two estimates into (6.19), we obtain that (6.19) is 
bounded above by 



(6.20) 



—a 00 

^— J^P (r > k \V Q = v) E [V? \V Q 



Ml 



k=l 
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Now repeat the same argument, but applied to the first term on the right 
of (6.15). Essentially, this can be viewed as one of the terms in the previous 
sum, namely the term k = 0. In fact, the previous argument may be repeated 
without change to obtain 

(6.21) anP n (v, E) < P (r > \V Q = v) E [V? \V = M] . 

Substituting these last two equations into the right-hand side of (6.15) and 
substituting (6.16) into the left-hand side of (6.15) yields 

(6.22) E [N U \V = v] = ^— J2^{r>k\V = v)B[V 1 a \V = M] 



< ^— { sup E w [ T ] E m [Vf] } , 

1 ~ P [we[-M,M] J 

which is (6.14). 

Finally observe that the quantity on the right-hand side of (6.22) is fi- 
nite. To this end, it is sufficient to verify that the set 93T is regular. For this 
purpose, observe that since {i] n } is an i.i.d. sequence, it follows by a slight 
modification of Lemma 5.1 (iii) that Wl is petite. Moreover, letting k denote 
the first return time of {V^} to Cm, then by Lemma 5.1 (i) and Theorem 
15.0.1 of Meyn and Tweedie (1993), we have that sup v ^ M E v [t K ] < oo for 
some t > 1. Then by definition, it follows that this last equation also holds 
for the augmented chain with 9Jt := Cm x {0,1} in place of Cm- Conse- 
quently, the conditions of Meyn and Tweedie (1993), Theorem 11.3.14 (i), 
are fulfilled and hence 9Jt is regular. 

Step 3: Finally observe that by summing the expectations studied in Steps 
1 and 2, we immediately obtain (6.5). □ 

PROOF of Theorem 4.1. By Proposition 6.1, it is sufficient to show that 
for some t G (0, 1), 



p 

-a [ 



(6.23) E 



"K-(u«)-l 

^ 1 {Vn>u} 



n=0 



V 

u 



U(\ogv), 



where U(z) := X^neN P*a(~ °°' z ) anc ^ PA is the marginal distribution of 
— log A. As explained in Section 4, U may be viewed under a minor conti- 
nuity condition as the renewal function of — S n = — X^=i 1°S^«> while the 
expectation on the left-hand side may be viewed as a truncated renewal 
function of the nonlinear process {V^}. 
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In the remainder of the proof, we will suppress the conditioning in (6.23). 
To prove (6.23), we first establish an upper bound and then a correspond- 
ing lower bound. For the upper bound, begin by observing (cf. (5.1)) that 



(6.24) 



\V n \ . . j_ (A n \D n \ + \B n \) 



\V n -l\ ~ \Vn-l\ 

Note by definition that |Vfc| > u* for all k < K(u l ). Hence 

(6.25) log^JMJ <log^ n + «-*^ n |Z? n | + |B n |)^, aJl n < 



Now introduce the random walk 



(6.26) :=J2 X i > where X i ■= lo § ( A i + ^(AnlDnl + \B n \)) 

i=i 



and S, 



(«) 



0. It follows from the previous two equations that log|V^| 



log | Vo| < S^K Since (Vq/u) = v, it then follows that 
(6.27) log \V n \ -\ogu< + log v. 

Now suppose that \V n \ > u. Then the left-hand side of (6.27) is positive, 

(u) 

and thus Sn > — log v. Consequently, 



(6.28) E 



'K"(tt*)-1 

E 1 {|V„|>«} 
n=0 



< E 



.n=0 



{Si u) >-\ 0g v} 



:= l/( u )(logv). 



To establish an upper bound, it remains to show that 
(6.29) lim sup U {u) (log v) < U{log v). 

To this end, observe by (H2) that E[exp jeX^}] < 00 for some e > (since 
we take the logarithm of the random variable A n \ D n \ + \ B n | in the definition 



u. 



of X\ u} ). Moreover, E[logA] < => E[A7" J ] < for sufficiently large 
Since can itself be viewed as a renewal function for the random walk 
—S^\ it follows that U^ u \logv) < 00. Since decreases monotonically 
to log Ai as u — > 00, it follows by a dominated convergence argument that 



(6.30) 



lim U iu) (logv) = E 



>— log V } 



n=0 



U(logv), 
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which establishes the required upper bound. 

Turning now to the lower bound, fix e > and choose N sufficiently large 
such that 



(6.31) 
Then 



E 



N 



^2 1 {S n >-\ogv} 



n=0 



— U(logv) 



< e. 





'K(u f )-1 




" N 


(6.32) E 


^ 1 {Vn>u} 


> E 


^ 1 {Vn>u} 




n=0 




_n=0 



NP (Kin*) < N) . 



To bound the first term on the right-hand side, note that 

V n > A n V n -i — \B n \, 
and by iterating this equation we obtain 

n— 1 n 

(6.33) V n > (Ai--- A n )V Yj II A j\ B *\> for a11 n - 

i=l j=i+l 

Thus, setting (Vq/u) = v, we obtain for any given n that 

(6.34) liminf — > [A\ ■ ■ ■ A n ) v := exp {S n + logv } a.s. 

Hence, if S n > — logv, then we will necessarily have liminf M _ >00 (F ri /u) > 1. 
Consequently, 



(6.35) liminf E 

u->oo 



N 



{V„>u} 



> E 



N 



n=l 



{5„>-log^} 



> U(logv) — e, 



where the last step was obtained by (6.31). Finally, observe by the definition 
of K{u % ) and (6.33) that for any s G (t, 1), 



(6.36) 



{K{u l ) < N) < P ( (Ai ■ ■ ■ A N ) V < n s ) 



(N-l N 
Y n >«"-'■ 
i=i j=i+i 

In the first term on the right-hand side, (Vq/u) = v, and so the probability in 
question reduces to P ((Ai • • • An) v < n s_1 ) = P (Sn + logv < (s — 1) log it), 
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and the latter probability tends to zero as u — > oo, as N is fixed and (s — 1) < 
0. Similar reasoning shows that the second term on the right of (6.36) also 
tends to zero as u — > oo. Thus we conclude lim^-^ P (K(u t ) < iV) = for 
any fixed N. Substituting this last equation and (6.35) into (6.32) yields the 
desired lower bound. □ 

PROOF of Theorem 4.2. First assume that (Vq/u) = v for some v > 1. 
Then it follows by (6.28) and Proposition 6.1 that 



Vo 

u 



E 



,n=0 



{V n >u} 



Vo 

u 



< U iu) {logv) + A(u), 



U {u) (logv) := E 



.n=0 



{5^>-log^} 



(6.37) E 

where 
(6.38) 

and is defined as in (6.26). In the definition of U^ u \ note that the terms 

inside the sum are obviously bounded above by 1 , „( u ) ^ , , and substi- 

{-SA ^<log-y}' 

tuting the latter quantity into the right-hand side of (6.38), we obtain the 
renewal function of the process {— S^}. Consequently, Lorden's inequality 
(Asmussen (2003), Proposition V.6.2) yields 



(6.39) 

where 
(6.40) 



+ 1 + 



a 



m 



m u := E[lW] and o\ := Var (X^) 



and is defined as in (6.26). 

Observe that these moments are actually finite. In particular, \x\ e dom- 
inates (log re) 2 for any e > 0. Hence, applying (H2) for a sufficiently small 
choice of e > 0, we obtain 



(6.41) 



E 



(!<»))' 



< const. • E 



< 00, for all u. 



Consequently, both the constants m u and a\ are finite. Substituting (6.39) 
into (6.37) yields 



(6.42) 



E 



N u 



Vo 

u 



< C\(u) log v + C 2 (u), 
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where (cf. (6.6), (6.40)) 



(6.43) 



Ci(it) :-- 



m,, 



and C 2 (u):= 1 + +A(u). 



Next recall that our primary objective is to study Ej> [iV u |&t u a(t-i)] • 
But by the strong Markov property, 

E S [iV u |ff T „A(T-i)] = Es, [iV u |y T J l {Tu<r} . 

Observe that in the dual measure, the process {V n } reverts back to its orig- 
inal measure after time T u , and if T u ft r then N u is zero. Thus, we may 
apply (6.42) to the right-hand side of the previous equation to obtain that 

(6.44) E 3 [JVuI^a^!)] < (Ci(iz)log(^ + C 2 («)) 1 { t u < t} , 

which is the required upper bound. 

Finally observe that as u — > oo, X^> decreases monotonically to log A 
a.s. and A (it) I 0, and hence as u — > oo, 



CAu) -> Ci := — and C 2 (u) -> C 2 := 1 + — 
m " ' m 



(6.45) 

where m := E [log A] and u 2 := Var (log A). 

PROOF of Theorem 4.3. By the strong Markov property 

Vn 

E TV,, -°= v :=H u (v), 



□ 



E 



V T 

1 u 

u 



u 



and by Lemma 4.1, 

(6.46) lim H u (v) = H(v) := UQogv). 

u— >oo 

Let v u denote the probability law of Vt u /u in the dual measure. Then our 
objective is to show that for any z > 1, 

rz rz 

(6.47) lim / H u {v)v-^dv u {v) = / H \v)v~* djl{v) . 

[Here we have taken the lower end point of the integral at one, since Vt u /u > 
1 for every u.] 

Next observe that 

(6.48) J H u {v)v^du u {v) = + R 2 (u), 
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where 



Ri(u) = J H u (v)v-td (v u (v) - ji(v)) , R 2 (u) = J H u (v)v-tdpL(v). 

Now by (6.42), (6.45), (6.46), and the dominated convergence theorem, it 
follows that as u — > oo, 

(6.49) R 2 (u) -^R 2 := J H(v)v~^djl(v). 

Finally, R\(u) — > as u — > oo using (6.42), (6.45), and the weak convergence 
of v u to p,. □ 



7. A generalization of Letac's Model E. While the results of the 
previous sections were obtained specifically for Letac's Model E, the methods 
are actually quite general, applying to rather arbitrary recursions subject to 
appropriate regularity conditions. However, if one departs from the setting 
of the previous sections, then one only obtains a recursive algorithm for 
computing the constant C, rather than the more elegant expression in terms 
of Letac's backward recursive equation. 

As an illustration of our method in a more general setting, consider the 
polynomial recursion 

(7.1) V = F Y {V) := B k V + B k _ l V { - k ~ 1),k + ■•• + B\V x l h + B , 

where (Bo, . . . ,B k ) is a sequence of positive- valued random variables and 
Y = (Bo, . . . , Bk-i). Note (7.1) may be written in the more suggestive form 

fc-i 

(7.2) V = AV + G Y (V), where A := B k and G Y (V) := BjV j/k . 

3=0 

Observe that when the coefficients (Bo, ■ ■ ■ , B k _{) are fixed, v~ 1 Gy(v) \. 
as v — > oo. Thus, for large V, the random function Gy(V) may be viewed 
as a remainder term, whereas the asymptotic behavior of the process is 
dominated in this case by the first term, AV. 

Starting with the SFPE (7.1), we can form the forward recursion 

(7.3) V n = A n V n -i + G Yn (V n -i), n = l,2,..., V = v, 

where {A n } is an i.i.d. sequence of random variables having the same prob- 
ability law as A, and {G Yn } is an independent sequence of random functions 
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having the same law as Gy and whose nth member will depend on A n . Es- 
sentially, (7.3) can be viewed as an extension of Example 3.2 of Section 3, 
where the term 11 B" in the SFPE (3.7) is now replaced with the random 
function Gy{v). With this modification, the entire discussion in Example 
3.2 adapts to this setting, yielding 



(7.4) 



C = 




As before Vq ~ v, where v is obtained from the minorization (Ai), while r is 
the first regeneration time and £ is the solution to the equation E[^] = 1. 
The expectation is computed in the ^-shifted measure, defined in the same 
manner as in (2.14), but with (A, Bo, . . . , -B^-i) in place of (A, B, D). While 
this new representation can no longer be viewed as backward iterations of 
the original SFPE, it nonetheless provides an algorithm with which the 
constant may be computed, and where the successive terms in the series 
become increasingly insignificant as n — > oo. 

To see that this series is convergent and to study the rate of convergence, 
using a slight variant of (5.34), it follows that the expectation in (7.4) is 
bounded above by 



n=l 



(7.5) ( E[V*\ + ]T E [(G Yn (V n ^)f] W P(r > n) 



where (I/O + (VO = (V0- 

To estimate E (Gy n (V n -i))^ as n — > oo, it is helpful to first observe 
that {V n } converges to a stationary distribution. Indeed, using a technique 
of Loynes (Goldie (1991), p. 161), it follows that for any positive constant c, 



'k-\ 



G Y (v) < c (/£ - 1)A max (l, -) + B , 



vi=l 



and thus Fy(v) := Av + Gy(v) satisfies 



k-i 



(7.6) F Y (v) < Ac + c(^V fc Bi max (l, ^) +B := A max (c, v)+B, 



i=i 



where B = Bq and A = Bj^j . The upper bound on the 

right-hand side is Letac's Model E, which itself converges to a stationary 
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distribution under the conditions of Goldie (1991), Proposition 6.1, namely 
under the assumption that £ > and that E [log(l V B)] < oo. This latter 
condition is subsumed in the first of the following hypotheses: 



(H 4 ) E [Bf ] < oo, % = 0, . . . , k - 1, for some f > £. 

fc-i 

! ! <() - 



K— 1 

(H 5 ) E[log(A+(^^)/2^) 



where B' is an independent copy of Bq. 

Thus we conclude that the Letac model given on the right-hand side of 

(7.6) converges to its stationary distribution. Next observe by (-H5) and the 
argument on p. 162 of Goldie (1991) that the backward iterates of the SFPE 
(7.1) converge and are independent of the initial value, and hence by Letac's 
principle (Lemma 2.1), the forward recursive sequence also converges a.s. to 
the same random variable, which we denote by V. Then Gy„(V^) converges 
a.s. to 

k-1 

(7.7) G Y (V):=Y,BjVV k , 

3=0 

where {Bj} is independent of V on the right-hand side. 

Finally, to obtain finiteness of the £'th moment of (7.7), note that the 
process {V^} is dominated from above by 

(7.8) V n :=A n V n - 1 + (A n c + B n ), n = l,2,..., V = V , 

which upon iteration yields % = Vq{Ai ■ ■ ■ A n ) Y™ =1 (AiC+Bi)(A i+1 ■ ■ ■ A n ). 
But given a < £, we can choose a constant c such that E[^4 Q ] < 1. Then a 
simple computation (using either Minkowskii's inequality in the case £' > 1 
or (3.4) otherwise) yields that E[t^f] is uniformly bounded in n for any a < 
£. Since {V n } dominates {V^}, it now follows, using in (7.7) the independence 
of V and {Bq, . . . , that 

k-1 

(7.9) lim E \(G Yn (Vn-i)f] =V E [ B f] E [K^l < 00 

3=0 

provided that £' has been chosen sufficiently small such that £'(&— l)/k < 
Then using the geometric recurrence of {V n }, it follows that (7.5) is finite; 
thus C is finite. 

For the above computations to make sense, it is critical that we can 
also verify the regularity properties of the Markov chain {V^}. But ip- 
irreducibility and the minorization and drift conditions all follow by very 
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minor modifications of the proofs given for Letac's Model E. Furthermore, 
transience follows in the ^-shifted measure, since under our assumptions 
V\ > Bq > 0, and then V n > Bq(A2 ■ ■ ■ A n ) f oo as n — > oo. Thus we have 
established the following: 

Theorem 7.1. Assume the polynomial model (7.1), and suppose that 
(H±), (H4), and (H5) are satisfied. Then 

(7.10) lim u f P (V > u) = C, 

u— >oo 

where C is a finite positive constant given by (7.4). 

It is also straightforward to obtain an upper bound along the lines of The- 
orem 2.2. The methods outlined above can be modified to analyze a variety 
of recursions of the general form V = AV + Gy(V) for an asymptotically 
"dominant" term AV and "remainder" term Gy(V), where v~ l Gy(v) I 
as v — > 00. For example, the composition / o h, where / corresponds to the 
SFPE in (7.1) and h is as in Letac's model (2.7), can be handled in this 
manner. 

8. Markov-driven processes and other extensions. 

8.1. Markov-driven processes and the ruin problem with invesments. Fi- 
nally, we consider the case where the forward or backward iterations in 
Letac's Model E are governed by an underlying Markov chain. In this case, 
it is possible to first apply Lemma 2.2 to the underlying Markov chain to di- 
vide its excursions into independent blocks, and to observe that the resulting 
process itself has the form of Letac's Model E. Thus, the theory developed 
in this paper may still be applied. 

Several examples of this type have been considered in Collamore (2009), 
including Markovian generalizations of the ruin problem with investments, 
perpetuities, and the GARCH(1,1) financial time series model. 

For example, in the ruin problem with investments, the goal is to deter- 
mine = P (sup nGZ+ C n > u) , where (cf. (1.12)) 

(8.1) C n :=B 1 + A 1 B 2 + (A 1 ---A r ^ 1 )B n , 

and in a Markovian setting it is natural to assume, instead, that (A n , B n ) = 
h(X n ), where h is a deterministic function and {X n } is a Harris recurrent 
Markov chain. Now by Lemma 2.2, the chain {X n } may be subdivided into 
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independent blocks of the form {Xx i _ 1 , ■ ■ ■ ,-Xj^-i}. Consequently, upon 
setting J* = J*(n) : = inf{i : i<Q > n} — 1, we obtain 

(8.2) £ n = Bo + + • • • + (i • • • A/*_i)5j* + (i • • • 
where 

M = A?q_i • • • -A-Ki-i, 

Bi = B Ki _ L + A Ki _- L B K% _ 1+ i H h (A Ki ^ ■ ■ ■ A Ki - 2 )B Ki -i, 

and lZ n = Bkj* + Akj* Bk^+i + • • • + {Arj* ■ ■ ■ A n -i)B n 

for all i = 0,1,..., where K-\ = 1. Now by a short argument, given in 
Collamore (2009), it can be shown that if we start the chain {X n } at a 
regeneration time, then it follows from (8.2) that C := sup ngZ+ C n satisfies 
an SFPE, namely, 

(8.3) C = max (M, AC) + B, 

where (A,B,M) = (A^B^Mi) and 

Mi := sup{Bj fi _ 1 + A jr ._ 1 B jr ._ 1+ i + ■ ■ ■ 

+ (A^.! • • • Aj-^Bj : < j < Ki} - B h for all i. 

Finally, it can be verified that — even without starting the process at a re- 
generation time — the asymptotics of P (C > u) obtained via the SFPE (8.3) 
are nonetheless the correct. For details, see Collamore (2009), Lemma 4.2. 

The main conclusion to be drawn from the above discussion is that the 
Markov-driven process also satisfies an SFPE having the form of Letac's 
Model E. However, if the Markov chain takes values in a general state space, 
as opposed to a finite state space, then it is not immediately clear whether 
the required moment conditions in (Hi) and (H 2 ) of this paper will actually 
be satisfied. For the ruin problem with investments, this issue has been ad- 
dressed in detail in Collamore (2009), Theorem 4.2 and Section 6.1. Namely 
set S n = ^r=il°§^> and consider the Gartner- Ellis limit from large devi- 
ation theory, 

A(a) := limsup - logE [e aSn ] , Va, 

and assume that the equation A(£) = has a unique positive solution. Next, 
set S n = Yli=i ^°sAi + log \B n \ and assume that both A(a) and 

A(a) := lim sup — log E 

n— >oo Tl 
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are finite in a neighborhood of £. Moreover, assume that {X n } satisfies the 
weak regularity condition (H3) and minorization condition (9Jt) given on 
pp. 1411-2 of Collamore (2009) (taking h(x) = e €X in (9Jt), for any e > 0). 
Together, these form one set of possible conditions suitable to apply the main 
results of Collamore (2009). Then applying our current results to identify 
the constants of decay, we arrive at the following extension of Theorem 2.1 
of Collamore (2009). 

Theorem 8.1. Assume the conditions of Collamore (2009), and assume 
that (A,B,M) satisfies the condition (H3) of this paper. Then 

(8.4) ^(u) ~ Cu~^ as u -)• 00, 

where C is given by (3.3) but with {(A4, Bi)}i £ i + replaced with {(Ai, Bi)}i £ z + , 
defined following (8.2), and where Bq is independent with Bq ~ £(B + ). 

8.2. Markov-driven Letac models. The previous discussion can be ex- 
tended by observing that (8.1) is the backward recursive sequence of the 
SFPE f(v) = Av + B, and we could equally well study the supremum of 
a backward sequence corresponding to a general Letac model, viz. f(v) = 
Amax(D,v) + B. As in the previous example, the limiting perpetuity se- 
quence for this more general model can be subdivided into i.i.d. blocks. 
Namely, letting Z n denote the nth backward iterate of the process, and as- 
suming that the driving sequence for these backward iterates is of the form 
(A n ,B n ,D n ) = h(X n ) for a Harris recurrent Markov chain {X n }, then just 
prior to regeneration, we have 



(8.5) 



V Ki -i = max 



(D^AiVK^-i) +B 



where 



Ki-l j-1 



Ai = A Ki 



1 ---A Ki _ 1 , Bi= Bj Q A k , 



j=Ki-i k=K i - 1 




Setting Z = lim^oo Z^ and taking the limit as i — > 00 yields the SFPE 



(8.6) 




max 



(D,AZ)+B, 



and under conditions analogous to those of Theorem 8.1, one naturally ex- 
pects this limiting tail distribution to correspond with that of limn^oo Z n . 
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We note that a slightly different SFPE will be obtained if, instead, one con- 
siders the supremum of the process {Z n } (involving a slight modification 
of the random variable M appearing in (8.3)). However, in either case, the 
resulting SFPE has the same form as Letac's Model E, and hence the theory 
developed in this paper applies. 

Forward iteration can be handled similarly, although the random variables 
(A,B,D) will be different than those obtained from backward iteration. 
Namely one obtains 

(8.7) V Ki -i = max (A, AVk^-i) + B t , 

where 

Ki-l K t -l 

A = A K% _ X ■ ■ ■ A Ki _!, Bi = ^ fij [J A k , 

i=Ki-i k=j+l 
Ki-l ( Ki-l k-1 Ki-l 

a = v { d * n a j - e b > n a i \ ■ 

k=Ki-! { j=k j=Ki-i l=j+l ) 

As a special case, forward iteration corresponds to a GARCH(1,1) process 
with Markov modulated innovations; see Theorem 2.2 of Collamore (2009), 
where the corresponding moment and regularity conditions are studied in 
a general Markovian setting for such processes. Again, Theorem 2.1 of this 
paper may be applied, although in all of these cases, it would be of interest 
to determine whether the pair (A, B) or triplet (A, B, D) appearing in the 
expression for the constant C can itself be characterized probabilistically in 
a more transparent way. 

9. Concluding remarks. In this paper we have introduced techniques 
that aid in the characterization of the constants C and £ for solutions to 
SFPEs in the stationary case. Using the results of this paper, we have con- 
currently developed statistical methods for joint asymptotic inference con- 
cerning C and £. This allows for the development of inferential methods for 
such quantities as Value at Risk and Expected Shortfall for GARCH (1,1) 
and related financial time series models, which are used extensively in ap- 
plied econometrics. Extensions of our work to continuous time processes and 
general Markovian processes are presently under study. 
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